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PREFACE 
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* 


1. OBJECTIVE 

The purpose of this report is to summarize the performance of work 
for the period 1 October 1982 through 30 March 1983, in compliance with 
Modification 19 to Article XXI of Contract NAS5-20682, entitled "Plasma 
Wave Experiment for ISEE-C (Heliocentric) Mission" dated 20 November 1974. 

The objective of this contract is to provide analysis of data from 
a scientific instrument designed to study solar wind and plasma wave phenomena 
on the ISEE-3 Mission. 

2. SCOPE OF WORK 

Project activities during this past six months have included successful 
return of data from the instrument, continuing analysis of all data, publication 
of results, and deposit in National Space Science Data Center of the data. 

3. CONCLUSIONS 

Not applicable. 

4. SUMMARY OF RECOMMENDATIONS 


Not applicable. 


1.0 INTRODUCTION 


The purpose of this report is to summarize the various activities 
and tasks accomplished on the data analysis phase of the contract during 
the last six months. 

2.0 WORK ACTIVITIES FOR THE SIX-MONTH PERIOD 


2.1 Research 

During this interval, there was much research activity involving 
the I3EE-3 Plasma Wave Investigation and the (unfunded) counterparts on 
ISEE-1 and -2. The paper entitled "Computer-Constructed Imagery of Distant 
Plasma Interaction Boundaries" by Greenstadt et , was published in 
Advances in Space Rearch (^, 7, 163, 1983). Several new papers were 
completed: "Plasma Boundaries and Shocks" (Russell and Greenstadt), 

submitted to Revs, of Geophys. and Space Phys ., February 1983; "Transfer 
of Pulsation-Related Wave Activity Across the Magnetopause: Observations 

of Corresponding Spectra by ISEE-1 and ISEE-2" (Greenstadt et al .) , 
prepared for the special issue of Geophysical Research Letters covering the 
AGU Chapman Conference on Waves in Magnetospheric Plasmas (Kona Coast, 
Hawaii, February 1983) and submitted to the journal, March 1983; "Science 
Return from ISEE-3 at Comet Giacobini-Zinner" (Scarf ^^.), to be 
published in the proceedings of the International Conference on Cometary 
Exploration (Budapest, Hungary, November 1982); and "The Interplanetary 
Shock Event of November 11/12 1978 — A Comprehensive Test of Acceleration 
Theory" (K.-P. Wenzel et al.), submitted to the International Cosmic Ray 


Conference {to be held in Bangalore, India, August 1983). Two other papers 
are in the process of being revised and readied for submission to the 
journals: "Collisionless Snock Lengths" (Mellott and Greenstadt), for 

Physical Review Letters ; and "Plasma and Energetic Particle Structure of a 
Collisionless Quasi-Parallel Shock" (Kennel ^^.)i prepared for Journal of 
Geophysical Research ." 

2.2 Other Activities 

During this six-month period. Dr. Scarf participated in the International 
Conference on Cometary Explorc:tion (held in Budapest, Hungary, November 1982), 
and the attached ISEE-3 manuscript was submitted for publication in the proceed- 
ings of the conference. Dr. Scarf also discussed the ISEE-3 mission to Giacobini- 
Zinner at the Inter-Agency Consultative Group meeting and at the meeting of 
the lACG Plasma Science Working Group. Dr. Scarf attended the ISEE Science 
Working Team meeting at Goddard Space Flight Center on Mar-rh 7-8, 1983, and he 
submitted ISEE-3 plasma wave data covering the first one and one-half crossings 
of the distant geomagnetic tail (October 1, 1982 through January 8, 1983) 
to NSSDC at this meeting. Mr. Greenstadt attended the A6U Fall Meeting held 
in San Francisco (December 1982) and presented a talk entitled "A Storm-Time, 

Pc 5 Event Observed in the Outer Magnetosphere by ISEE-1 and -2: Wave 
Properties". In February he traveled to Kona Coast, Hawaii, to participate 
in the Chapman Conference and to present a talk on "Transfer of Pulsation- 
Related Wave Activity Across the Magnetopause: Observations of Favorable 
Conditions by ISEE-1 and -2". 
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INTRODUCTION , “ 

The quadrenniun has seen an explosive growth 
in our knowledge and understanding of the various 
plasma and magnetic field boundaries in the 
terrestrial magnetosi^ere. Nowhere is that more 
evident than at the magnetopause and bow shock. 
In the 1975 and 1979 quadrennial reports the 
magnetopause was covered in a single paragraph 
each year. The bow shock received one paragraph 
in 1975 and none in 1979! The reason for the 
resurgence of interest in these boundaries was 
the availability of new and exciting measurements 
ffom the ISEE-1 and -2 spacecraft. Not only did 
these spacecraft carry sensitive high-time 
resolution, three-dimensional plasna 
instrunentation as well as high-time resolution 
and accurate magnetic and electric field data, 
but also the variable separation of the two 
spacecraft allowed the velocity of structures to 
be measured and thereby allowed time profiles to 
be converted to spatial profiles. Simultaneously 
mnerical simulationists were benefiting from 
larger and faster computers and to development of 
ever increasingly sophisticated codes. During 
this period the simulation field began to 
blossom. The ccxnbination of good data and 
realistic models led to not only empirical 
knowledge but theoretical understanding of many 
of the processes at work. 

The quadrenniun interval 1979 through 1982 
spanned the transition from exploratory work with 
preliminary data from ISEE-1 and 2 [e.g. Bame et 
al., 1979; Russell and Greenstadt, 1979; 1981] to 
detailed investigations with fully calibrated and 
updated measurements fVom these satellites and 
ISEE-3 as well. Published ccmtributions 
therefore included the last of the ISEE 
preliminary papers, annomcements of new 
discoveries, and follow-up accounts of plaana 
booidary phenomenology,— plus reports based on 
Voyager and Pioneer-Venus data. 

Progress in the field was aided by many 
workshops and conferences devoted to the 
magnetospheric boundaries. The first of these, 
the Chapnan Conference on Magnetospheric Boundary 
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Layers, held in Alpach Austria in June 1979 
covered mainly the magnetopause and boundary 
layer. The proceedings of this conference has 
been published as an ESA special publication, 
There were two coordinated data analysis 
workshops, CDAW-3 and CDAW-i|, devoted to the bow 
shock and the magnetopause as observed by ISEE-1 
and -2. These have been described by Qgilvie 
[1982] and Paschmann [1982], respectively. ‘The 
latter workshop spawned a series of papers in a 
special issue of the Journal of Geophysical 
Research [April, 1982]. A Gordon Conference on 
collisionless stocks was held in June 1981 but as 
is their policy no proceedings were published; 
however, a workshop on upstrean particles and 
waves held at JPL produced a special issue of the 
Journal of Geophysical Research: ISEE: Upstream 

Waves and Particles, Vol. 86, No. A6, June 1, 
1981, which contains many of the articles cited 
here, including overviews by Tsurutani and 
Rodriguez [1981] and Kennel [1981]. 

In this review we will proceed fl*om the 
outside in. We will discuss first the bow stock 
and foreshock and say a little about 

interplanetary stocks. Then v#e discuss the 
magnetosheath, magnetopause and bouixlary layer. 
After a section on reconnection we treat the 
plasma and neutral sheets, polar cusp and the 
injv-ction of plasma into the inner magnetosphere. 

SHOCKS 

General 


Broadly speaking, studies have divided 
themselves into two categories: careful 
documentation of more or less anticipated 
features of the stock using the high-quality 
instrumentation and two-point measurements of the 
lSEE-1,2 spacecraft, and intensive investigation 
of the [quasi-parallel] foreshock inspired by the 
discovery of the array of return ion 
distributions reported in the last lUGG interval 
[Gosling et al., 1978]. Both cat^ories have 
benefited from reawakened interest in, and 
support of, theoretical calculations and 
large-scale computer simulations of stock and 
foreshock phenomena. As might be imagined, 
routine docunentation is preceding slowly because 
of the care required in quantifying stock 
features reliably and because of f1*equent 
diversion of a limited nunber of researchers to 
the exciting area of foreshock dynanics. 

Figure 1 illustrates conceptually the items 
of stock structure that have been the principal 
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foci of activity to the present. Figure 1 is an 
adaptation of a sketch used in an earlier review 
prepared at the beginning of this lUGG interval; 
the hatched areas indicate these combinations of 
shock geometry and strtjcture that have been the 
subjects of intense and detailed Investigation. 
We see that substantial areas of the figure are 
uishaded, and that only a minor fraction of the 
shock’s features have been studied. 
Nevertheless, much has been learned and a 
sii)stantial foundation for continued study has 
been laid. We note here that an additional area 
of investigation has been opened by 
identification of electrostatic shocks in the 
magnetosphere between 2.5 and 7 Rc on auroral 
field lines [Mozer, 1981]. ^ 

Hydromagnetic Modeling 

Refinement of bow shock location and shape 
continued with efforts to compare messuranents 
and models at the terrestrial planets [Slavin and 
Holzer, 1981; Slavin et al., 1983; Mihalov et al. 
1982], and the analytic description of the bow 
shock system was advanced by Zhuang and Russell 
[1981] in a comprehensive study of magnetosheath 
thickness. Encounters with the Jovian bow shock 
were described by Lepping et al. [1981]. Harvey 
et al. [1981] and Bonifazi et al. [1980b, c, 
1982] foind they could routinely use simple MHD 
formulas for approximating the bow shock location 
needed in their studies of upstream conditions. 

Bow Shock Front or Thermalization Layer. 

By either term above we mean that feature of the 
shock wherein a substantial portion of the solar 
streaning energy of any plasma component is 
converted into thermal, although not necessarily 
isotropic or roaxwellian, energy. The thicknesses 
and the diagnostics by which the layer is 
defined, as well as the kinetic details of the 
conversion are the subjects of inquiry. 

An overview of the "typical” 
quasi- perpendicular shock was given by Greenstadt 
et al. [1980a], in which the relationships of the 
various diagnostics to each other could be seen 
at medium resolution, by current standards. Rapid 
heating of electrons in the magnetic foot and 
ramp of the shock was apparent, along with 
generation of a secondary distribution of 
accelerated protons, vAiile thermalization of both 
proton distributions was seen to occur through a 
more prolonged series of magnetic oscillations 
behind the principal shock ranp and overshoot. 
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Ihe overshoots themselves were studied by Russell 
et al. [1982a] and Livesy et al. [1982], who 
fouid that the overshoots grow with Mach nunber 
above critical [M>3]; this result is consistent 
with that of simulations showing ^ association 
of overshoots with trapped, reflected particles 
constituting the secondary distribution of 
protons present in the shock fy*ont above M'S 
[Leroy et al.,1981; 1982], 

High resolution views of electron heating at 
the bow shock were described by Bane et al. 
[19'^9]* Further details of electron distributions 
typical of the q-perpendicular shock have been 
displayed by Feldman et al. [1982]: Electron 

distributions are skewed and anisotropized by the 
shock potential, producing distribution envelopes 
with "bunps", i.e. f>"ee-energy components defined 
by nonmonotonic df/dE, offset along the anbient 
field toward or away from the shock, depending on 
location up- or downstream. I^ch distributions 
can be expected to excite the plasma 

insta ilities responsible for therroalizing the 
particles, as diseased in the Feldman et al. 
paper and el£^rated theoretically by [Thomsen et 
al., 1982]. A careful analysis of the energy gain 
of electrons in the observed reference f1"ame has 
been undertaken by Scudder and Gocxlrich [1982] , 
and an extensive treatment of microinstabilities 
in the bow shock has been prepared by Wu et al. 
[ 1982 ]. 

High-resolution observations of ion heating 
and related behavior in the bow shock emerged 
with a paper by Paschmann et al. [19821, who 
displayed tw^imensional ion velocity 
distributions through the foot, ramp, and 
overshoot of a nearly perpendicular shock, at 
6-sec resolution. Figure 2 [taken from Figure 1 
of Paschmann et al.] illustrates schematically 
the ion behavior inferred from the spacecraft 
data, which agrees well with the results 11*0111 
simulations already cited [Leroy et al., 1981; 
1982]. The first simulation report [Leroy et 
al., 1981] showed that the "overall shock 
structure consists of several distinct regions 
whose properties are closely cc^r^ected to the 
dynamics of the reflected ions." Ihe second 
simulation report [Leroy et al., 1982] described 
the effects on perpendicular shock configuration 
of a wide range of parameters: Scaling, 

stationarity , and reflection were all studied as 
functions of beta and Mach nunber. Perhaps the 
most striking result was the variation in numbers 
of ions reflectioned, ft-om none at subcritical 
Mach nunbers to 30-^oi at M(AIfven] )12-13. 
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A related look at observational details of 
ion dissipation with the ISEE spacecraft was 
provided by Formisano and Torbert t1982]» who 
correlated large>«nplitude electrostatic waves in 
the frequency range 2-512 Hz with the presence of 
doii>ly-peaked ion energy shock profiles. Ihe 
results l«id support to the long-standing belief 
that ion-ion streaming produces the instability 
generating ion acoustic waves responsible for. ion 
heating in supercritical collisionless shocks. 

Ihe general picture of quasi-perpendicular 
shock structure has been enhanced by a 
description of the variable plasma wave 
turbulence at different planetary shocks, where a 
progressively more distinct spectral peak between 
the Bmeman and ion plasna frequencies was found 
in going from Venus to Saturn [Scarf et 
al.,198l]. It was speculated that this 
progression followed the increase in Hach number 
with distance in the solar wind away from the 
sun. At an opposite extreme, Russell et al. 
[1982b] have been attempting to analyze Earth's 
bow shock in its simplest form at low Mach nonber 
and low beta. They find that laminar, 
quasi-perpendicular shrcks defined by these 
conditiOTS have thicknesses close to the ion 
inertial length c/upi and that thickness 
increases as 6 Bn drops toward 55 deg. Figure 3 
sketches the rough dependence of Q-perpendicular 
shocks on Mach nimber and ar<;le. 

Study of quasi-parallel ^hock structure is 
progressing very slowly, A great amount of data 
has been examined, but almost nothing has been 
published that doesn't center on the associated 
foreshock, reviewed separately below. Recent 
examination of foreshock ions, however, has found 
evidence of specularly reflected ions in or just 
ahead of quasi-parallel transition, implying the 
presence in q-parallel structures of potential 
layers capable of causing such reflection 
[Gosling et al., 1982]. Fresh simulations of 
q-parallel shocks [Quest et al., 1982] have 
generated specularly reflected ions and have also 
found strong ion heating, stronger thai electror 
heating, at low HB;'h nunber, for a specific 
parameter set chosen to be similar to that of the 
cases observed by Gosling et al., [1982]. Ihus, 
simulations show considerable promise of 
improving ou* understanding of quasi-parallel 
structures in the near future, with return ions 
becoming valuable as diagnostics of q-parallel 
structure. Ihe first quantitative attack on 
Q-parallel macro- structure with lSEE-1, 2 was 
reported by Greenstadt et al. [1982a], who 
computed the correlation between large- amplitude 
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Di^netic {MJlsations at the two spacecraft at the 
outer edge of the shock structure, finding a 
sharp drop from high to low correlation at a 
satellite-separation of about 1000 Ion. The drop 
appeared to be compatible with Larmor radii of 
protons typically observed in Uie foreshock just 
upstream from the outermost pulsations. 

Figure summarizes, by shading portions of 
a familiar table of structural designations , ‘the 
types of shock structure where the principal 
activity has been taking place. 

Foreshock 


Particles . Starting with the well-defined, 
energy-dependent electron boundaries found near 
the field line tangent to the bow shock [Anderson 
et al.,19793, foreshock, i.e. return, particles 
have ccimanded major attention. The 
distinguishable distributions of 1-40kev return 
ions discovered earlier — beams, diffuse, and 
intermediate [Gosling et al. ,19793— have been 
tied by simultaneous IMF directions to 
Q-perpendicular, Q-parallel, and transition shock 
structures, in the same order [Greenstadt et 
al., 1980b]. The characteristic of the reflected 
beans appeared tc be consistent with a model 
proposed earlier by Sonnerup [19693, based on 
preservation of magnetic moment in the reflection 
];rocess [Paschmann et al., 1980] ar. all the 
various particle classes have been characterized 
in phase space [Eastman et al., 1981; Sentman et 
al. 1981a; Gurgiolo et al., 1981; Formisano et 
al., 1980], and the deceleration of the solar 
wind by diffuse distributions has been noted 
[Bame et al., 1980] and confirmed by comparison 
of data from two satellites [Bonifazi et al., 
1980c]. At last, theoretical instability 
computations have been performed using realistic 
ion and electron distributions in the solar wind 
[Gary, 1981; Gary et al.,1981; Sentman et al., 
1981b; Feldman et al., 1982b], and consequences 
of these are being tested with further data. 

Additional interest has centered on higher 
energy return protons in the 30kev to IMev range, 
particularly on their detection and occurrence 
far upstrean at ISEE-3 as well as near the shock 
at ISEE-1 [Anderson, 1981; Clceckler, 1979; 
Ipavich et al., 1979a,b; 1981a, b; Sanderson et 
al., 19813. Alrhoi^ some sunward streoDli^ ions 
have been attributed to the magnetosphere 
[Scholer and Hovestadt, 19813, anisotropies and 
spectr.l shapes of most of these ions, together 
with their link to the quasi- parallel bow shock 
via correlation with the ap^opriate IMF 
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directions [Scholer et al., 1979a, b; 1980a, b; 
1981a] have awakened interest in seeking 
observational evidence near the earth's shock 
that cosmic rays can be fashioned out of the 
plasma background by reflection, or at least 
ejection, fl^am relatively weak shocks like the 
earth's. Association of 40kev and 30Mev protons 
and, indeed, the entire range of ion energies 
from IKev to 1.6 Hev, with both upstream -axl 
dowistream regions of interplanetary ' shocks 
[Gosling et al., 1980-1981; Even^^n et al., 
1982], most of which are even weaker than the 
earth's at 1AU, have strengthened this interest, 
to be discussed again in a later paragraph. 
Alpha particles have also been found to 

participate in the reflection of solar wind 
particles from the bow shock [Scholer et al., 
1981b]. 

Waves . Athough foreshock waves are derived 
from interactions of foreshock particles with the 
solar wind, they were the phenomena from which 
the foreshock was originally defined, because of 
their ease of detection and processing compared 
to backstreaning particles. Waves have continued 
to be the pioneer tool in investigating upstream 
effects. In particular, they have served to 
demonstrate the existence of foreshocks of other 
planets. H^ppe and Russell [1981] showed the 
universality of the various foreshock ULF waves 
at Mercury, Venus, and Jupiter, including 
discrete wave packets and the incipient 

shock-like wave gradients underlying them at 
Venus. Indeed, the dependence of upstream wave 
periods on IMF magnitude at Earth [Russell and 
Hoppe, 1981] was extended to a coonon 

relationship for the waves at four planets, 
implying resonance with beams of ions of 
essentially the sane energy at each of the 
planets, and perhaps in other astrophysical 
systems as well [Hoppe and Russell, 1982a]. The 
wave foreshock concept was also extended to 
interplanetary [IP] shocks, whe*‘e VLF 

electromagnetic and f^asma waves were fomd to 
precede IP shocks for many hours ahead of actual 
shock arrival ii the local geometry was 

Q-parallel with respect to the IMF overtaken by 
the front [Kennel et al., 1982; Greenstadt et 
al., 1982a]. 

Wave properties in the earth's foreshock 
were also explored further. Wave packets 
attached to the largest amplitude ULF waves were 
found to be vriiistlers by correlating ISEE-1 and 
-2 data [Hoppe and Russell, 1980], and small 
waves in the foot of ^ Q-perpendicuLar shock, a 
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"mini foreshock”, so to speak, were shown to be 
consistent with phase and group velocity 
properties of v^istlers in the local solar wind 
[Greenstadt et al., 1931]. Upstream plaana waves 
were described by Filbert and Kellogg [1979] and 
Gurnett et al. [1979]. 

Particles and waves . The physics of the 
foreshock lies, of course, in neither particles 
nor waves, but in the relationships between them. 
FiTidamental to these relationships is the 
correspondence between particle and wave types 
[Paschraann et al., 1979; Hoppe et al., 1981; 
1982a,b], according to which beans, intermediate, 
and diffuse ion distributions are associated 
respectively with small, 1Hz whistlers, larger 
transverse narrow band, .1<-.01Hz waves, and still 
larger [6B'B/2], compressional , .1-.01Hz waves 
and connected whistler-wave packets. Figure 5, 
from Hoppe and Russell [1982b] sunmarizes this 
correspondence. Dependences of properties of the 
ion distributions, e.g. d^sities, velocities, 
temperatures , on position within the foreshock 
and its wave structr’-e have been extensively 
developed by Bonifazi and Moreno [1981a, b], and 
Bonifazi et al., [1980a, b, 1982]. Instabilities 
caused by beams [Gary et al., 1981] and diffuse 
distributions [Sentman et al., 1981b] have been 
proposed to explain the associated waves, as 
mentioned earlier, but neither theory has 
correctly provided all the observed properties of 
the waves. In any case, skewed foreshock electron 
distributions offer a better and more likely 
explanation for the 1Hz waves than the ion beams 
[Feldman et al., 1982a], which are not always 
seen with these waves [Hoppe et al., 1982]. In 
the VLF range, narrow- banded ion acoustic waves 
are associated with backstreaming protons under 
about 1.5keV, while electron plaana oscillations, 
whistlers, ion acoustic waves, and low frequency 
electrostatic waves accompany electrwis between 
.2 and I.SkeV [Anderson et al., 1981]. The 
upstream plaana waves are closely allied with 
spatial gradients of foreshock particles and 
bursts of wave noise appear at the edges of 
particle enhancements [Anderson et al.,op.cit.; 
Parks et al. , 1981]. 

Special attention has focused on docunenting 
the properties of foreshock ions in hooe of 
discovering the origin of the diffuse ion 
distributions and ULF waves accompanying them. 
The puzzle of the foreshock centers on the source 
and maintenance of the ULF/diffuse foreshock, 
which is doii)tless the product of some particle 
interaction, since the waves cannot by themselves 
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propagate upstream . It has been suggested [ 1 ] 

that the diffuse ions evolve out of an 

interaction of the reflected beams with the solar 
wind that produces waves which in turn disrupt 
the beams and scatter the beam particles into 
their diffuse distributions [Baroe et al., 1980; 
Gary et al., 1981; Pasclinann et al., 1981]| and 
[2] that the diffuse ions derive from Inherent 
properties of the q>parallel shock structure with 
whose geometry their occurrence is so intimately 
connected [Greenstadt et al., 1980b; Eastman et 
al., 1981]. Ihe first hypothesis is supported in 
part by weakening of the beams with distance 
[Bonifazi and Moreno, 1981b]; the second by the 
association of specularly, as opposed to magnetic 
moment-conserving, reflected ions specifically 
with q-parallel structures [Gosling et al., 1962; 
Quest et al., 1982]. It may reasonably be 
expected at this time that resolution of the 
problem will involve a combination of both 
models, facilitated by new theoretical approaches 
differentiating source populations [Schwartz et 
al., 1982]. 

Among the most exciting results in the 
foreshock has been the developnent of both theory 
and observation of the higher energy, ilOkeV-IMeV, 
protons. These have been treated by Monte Carlo 
simulation [Ellison, 1982] and as products of a 
Fermi process in which the upstream ULF waves 
serve as reflecting and scattering centers, 
together with the shock downstrean, to produce 
incipient cosmic rays [Terasawa, 1979; 1981]. A 
nunber of observations have supported the 
energization of return ions from the shock by 
repeated reflection between the scattering 
centers upstrean and the shock or other 
scattering centers downstream [Scholer et al., 
1979b; 1980a]. The upstream centers appear to 
scatter protons with a mean free path of about 4 
Rp and to lie within about 30 Rr of the bow shock 
[Scholer et al., 1980b]. A self-consistent model 
of ion energization using reflected beaas as 
feedstock to generate waves, which in turn 
reflect and scatter the ions until they leave a 
•'free escape” boundary beyond 10 or 20 Re 
upstream iVom the Q-parallel shock [Lee et al., 
1981; Lee, 1982], has succeeded in predicting 
energy spectra consistent with measured 
distributions [Ipavich et al., 1981a, b]. 

Interplanetary Shocks 

The growing theory and evidence that energetic 
ions are produced out of the thermal plasma by 
the earth’s bow shock has stimulated fresh 
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attention to IP shocks, where sup'athermal 
particles, both electrons [Potter, 198 u and ions 
[Gosling et al., I960, 1981] have been recorded, 
along with precursor plasma wave noise that might 
be indicative of forward particle escape [Kennel 
et al., 1982; Greenstadt et al, 1982b]. The 
structure of IP shocks has enjoyed renewed 
scrutiny [Pesses et al., 1979; 1981; Russell and 
Greenstadt, 1981], and an interesting instance of 
bidirectional electron streaning suggestive of 
closed field lines in the driver gas has been 
described [Bane et al., 1980]. Finally, coronal 
transient phenomena have been reviewed generally 
by Dryer [1982], and the interaction of IP shocks 
with the bow shock, magnetosheath, and 
magnetopause have been modelled and observed 
[Zhuang et al., 1981; Winterhalter et al., 1981]. 

Magnetosheath 

As usual , the magnetosheath has been the 
focus of massive inattention relative to other 
shock-related siAsjects, with a few notable 
exceptions. Early evaluation of icxi and electron 
measurements denonstrated the existence of 
’’quiet” and "disturbed" local states in the 
sheath, the latter distinguished by the presence 
of energetic ions, 3-40KeV, long period density 
fluctuations, and turbulent flow [Asbridge et 
al», 1978; Qgilvie and Scudder, 19793. The 
disturb^ state is probably the same as that 
identified with protons of energy > 100 KeV by 
West and Buck [19763. The average electron heat 
flux in the sheath was found to be double that in 
the solar vlnd and directed away from the shock 
[Ogilvie ans Scudder 19793. In a broad look at 
the magnetosheath, Crocker et al. [19813 compared 
the appearance of energetic ions in the sf.eath 
with IMF orientation to produce patterns of these 
ions downstrean. Because of the known correlation 
of the ions with enhanced magnetic sheath 
turbulence, ion patterns probably approximated 
the patterns of downstream Q-parallel structure, 
or so the comparisons suggest: IMF across the 
solar wind flow removed the ions to the outside 
lanks of the sheath; stream angle IMF placed them 
mostly in the morning sheath; IMF parallel to the 
solar wind filled the subsolar sheath with 
energetic ions. Figure 6, from the Crooker et al. 
paper, illustrates these cot figurations. The 
possibility that trapped, magneto spheric ions may 
provide some of the energetic ions found in the 
sheath was argued by Reiser et al. [1982], and a 
generation mechanisn for lion roars in the sheath 
was proposed by Thorne and Tsurutani [19813. 
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"Festoon”>shaped electric field emissions between 
0.1 and 4 KHz in f-t spectrograms have been 
id&itified as doppler-shifted ion acoustic waves 
convected through the sheath from the bow shock 
[Gallagher, 1982]. 

Magnetospheric Boundaries 

Magnetopause . Ihe anoint of research done on .the 
magnetopause over the last few years has been 
great and we cannot do it Justice in the space 
provided here. Readers interested in further 
details are referred to the various excellent 
reviews available. Fairfield [19793 simnarizes 
work on the magnetopause prior to the ISEE 
results. Paschraann [1982] provides an up-to-date 
review of the ISEE results and Sonnerup [19793 
reviews the theory of reconnection which is a key 
process at the magnetopause. Finally, Cowley 
[19823 provides a detailed and up-to-date review 
of reconnection at the magnetopause, flux 
transfer events, and the boundary layer in the 
guise of examining the causes of connection. 

Location . There is very little work being 
done on the average location and shape of the 
magnetopause at present. However, magnetopause 
motions are still of quite some interest. One of 
the puzzles about the shape of the magnetopause 
was why it appeared to be symmetric despite the 
fact the J X B forces would cause an additional 
deflection of The post shock solar wind. Zhuang 
et al. [1981b] developed an analytic model of the 
magnetosheath to examine this problem and found 
that, while the flow deflection did take place in 
front of the magnetosphere, the thermal and 
magnetic pressure of the magnetosheath maintained 
a nearly symmetric magnetosphere. 

The interaction of interplanetary shocks 
with the magnetopause was exanined by Grib et al. 
[19793 and Zhuang [1981a]. The observed response 
of the magnetopause to shocks was well explained 
by these models. Another cause of magnetopause 
motion was examined by Crooker and Siscoe [19793* 
They found that substorms caused large excursions 
of the tail magnetopause. 

One postulate that has attracted some 
vociferous support but little experimental 
evidence is that of impulsive injections through 
the magnetopause [ef. Heikkila, 1979> 1982a]. 
The advocates of this model suppose that over 
dense regions or blobs in the solar wind are 
blown against (and through) the magneto Tiuse. 
These models Ignore the fact that the stagnation 
streanline paints the entire magnetopause. The 
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si^solar magnetosheath spreads the normal 
stresses over the entire magnetopause. 

Gasdynanic simulations of the solar wind 
interaction with the magnetosphere have been 
available for many years, but only recently have 
there been magnetohydro dynamic models. Hcst 
recently these models have developed in three 
dimoisions [Leboeuf et al., 1981; Wu et al., 

1981] . Presently they do little more than 

determine the location of boundaries and the flow 
field aid field line geometry around the obstacle 
but even this is useful especially at low Mach 
nunbers where the gasdynanic solutions are 
expected to be inappropriate. 

Motion and Structure . Hie earliest ISEE 
measurements revealed the magnetopause to be in 
irregular and constant motion [Russell and 
ElphiCy 1978; Paschmann et al., 1978; Elphic and 
Russell, 19793. As shown in Figure 7 the 
velocities ranged from kilometers per second to 
hundreds of kilometers per second [Berchem and 
Russell, 1982a] but the thickness was much more 
constant at about ^00-1000 km or a few ion gyro 
radii. One of the surprises of the magnetopause 
was the variation in structure in short distances 
along the boundary. During a rapid 
interplanetary shock-induced motion of the 
magnetopause the currents on the magnetopause 
cannot be simply described in terms of MHD 
discontinuities. Although occasionally the 
magnetopause has the magnetic and plasma 
signatures of a rotational discontinuity 
[Paschmann el al., 1979» it also can be found 
with a clear tangential discontinuity signature 
even though the magnetosheath and magneto spheric 
fields are antiparallel (Papanastorakis et al., 

1982] . Contrary to coranon belief the magnetic 

field rotation is not controlled by the gyration 
of ions or electrons so that it follows a 
particular path dependent on the sign of the 
normal component of the field crossing the 
boundary. Rather the path the magnetic field 
follows from the magnetosheath to the 
magnetosphere is the shortest path [Berchem and 
Russell, 1982b]. Computer simulations have been 
performed of both the tangential discontinuity 
and rotational discontinuity by Lee and Kan 
[1979; 1092]. These originally favored the 

electron polarization but recent simulations 
[Swift and Lee, laipublished manuscript, 1982] 
reproduce these new observations. 

Hux Transfer Events . Frequently when 
ISEE-1 and are near the magnetopause, the 
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magnetic field will oscillate in a manner that 
resembles that occurring during a ms^netopause 
crossing but with some significant differences. 
Ihse differences suggest that a tube of magnetic 
flux in the magnetosphere has reconnected with 
some magnetosheath magnetic field and is being 
pulled tailward [Russell and Elphic, 1978, 1979; 
Ell^ic and Russell, 1979]* A sketch of the 
configuration of an FTE is shown in Figure 8. 
During these events the low energy ' plaana 
resembles the magnetosheath plasma and the 
energetic particles resemble the magnetospheric 
population. Daly et al. [1981] have shown that 
the energetic ions are streaning out of the 
mangetosphere. The electron signature is, 
however, more confusing. These results were 
confirmed by Scholer et al. [1982a]. Paschmann 
et al., [1982] have exanined the over-pressure in 
the flux transfer event and find that it is equal 
to the Maxwell stress imposed by the twist and 
draping of the magnetic field around the tube. 
Modeling the leakage of particles out of the 
magnetosphere through an FTE has shown that the 
ISEE observations can easily be replicated thus 
lending further support to the FTC interpretation 
[Speiser and Williams, 1982]. Initial] y all 
FTE's had the same signature in the component of 
the magnetic field normal to the magnetopause, 
outward then inward. However, Rijnbeek et al., 
[1982] have discovered reversed FTE's on August 
9, 1978 a period when Sonnerup et al., [1981] 
showed that ISEE was probably observing 
steady-state reconnection below the merging line. 
This is very important for two reasons. First, 
it shows that the sign of the FTC signature 
reverses from north to south across the merging 
line. Second, it shows that the magnetosphere 
may be undergoing steady-state reconnection and 
patchy reconnection simultaneously or nearly so. 
The statistical accuracy of these results have 
been extended using three years of data by 
Berchem and Russell [unpublished manuscript, 
1982] who show that the magnetospheric equator 
essentially divides the FTE signatures into two 
groups, normal and reversed. This pattern is 
consistent with FTE's being created by 

reconnection at the magnetospheric equator and 
the being pulled poleward away from the equator 
by field line tension and magnetosheath 

convection . 

Remote Sounding of the Magnetopause. 

Williams [1979, 1980] and Williams k al., [1979] 
note that the near presence of the magnetopause 
is associated with gyro phase asynmetries in the 
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3-D energetic ion observations. In other words, 
particles which intersect the magnetopause (or 
some boundary near the magnetopause) in their 
cyclotron motion about a field line appear to be 
lost. While Williams and coworkers originally 
identified this boundary with the magnetopause, 
it need not be so. It is perhaps more correct to 
call it the trapping boundary which could, for 
example, be the boindary between open and . closed 
field lines. This technique has also be used to 
measure magnetopause velocity [Fahnenstiel , 1981; 
Fritz and Fahnenstiel, 1982; Fritz et al., 1982] 
but because the technique, at least as applied to 
ISEE-1 data, returns a magnetopause location once 
every 36 seconds care must be exercised in the 
interpretation of the data. K>st recently Daly 
[1982] has critically exarained the remote 
sounding technique and concludes that the simple 
absorbing wall model overestimates the distance 
to the boundary. 

Boundary Layer . The botndary layer is 
perhaps tni outstanding enigma of the 
magnetosphere. It is a region of density and 
temperature intermediate between that of the 
raagnetosheath and that of the magnetosphere just 
inside the magnetopause. The thickness of the 
boundary layer is extremely variable [Eastman and 
Hones, 1979]. It also often is flowing away from 
the sun but occasionally flows the other way. 
Sonnerup [1979, 1980a] has developed a simple 
model of a viscous boundary layer coupled to the 
earth by field-aligned cur'*ents. Field-aligned 
currents certainly are present but they don't 
obey the simple Sonnerup model [Sckopke et al., 
1981; Hones et al., 1982]. In fact, the boundary 
layer studied by &kopke et al. [1981] seems to 
be associated with the magnetospheric roots of 
flux transfer events [Cowley, 1982]. This should 
not be surprising since Paschmann et al. [1982] 
find that flux transfer events in the 
magnetosheath contain a mixture of magnetosheath 
and magnetos;*ieric plaana having temperatures 
higher and densities lower than the 
magnetosheath. Such properties are very similar 
to those of the boundary layer. Reiff [1979] has 
reviewed the properties of vdiat is thought to be 
the low altitude extension of the bomdary layer. 
The low altitude measurements reveal two distinct 
plasma populations on the dayside: one 
predominatly on open field lines with low 
electron temperature (“15 eV), decreasing proton 
energy with increasing latitude and correlated 
with the IMF; the other on closed filed lines 
with a V-shap^ ion distribution and an electron 
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temperature of about 50 eV wanner than the above 
layer but cooler than the plaaroa sheet. 

Miyiellaneous . A very attractive postulate 
^ut the location of the merging line was put 
forth by Crooker [1979a,b; 1980] who assumed that 
merging was most likely to occur at t1x>se places 
on the magnetopause vAiere the magnetospheric and 
magnetosheath fields were exactly antiparallel. 
The resulting merging line depended bn the 
relative orientation of the magnetospheric and 
magnetosheath magnetic fields and was in the 
equator and passed through the sii>solar point 
only when the interplanetary magnetic field was 
exactly southward. Otherwise the merging line 
was at high north and south latitudes. However, 
observations of steady-state merging [Sonnerup et 
al., 19813 and FTTI's, as discussed above, support 
a near equatorial merging line at all tii^s. 

If reconnection is taking place, the plasma 
flows into the magnetopause from both the 
magnetosheath and the magnetosphere. This is 
equivalent to a tangential electric field. Such 
a tangential electric field has been reported by 
Mozer et al. [1978, 19793 and Fahleson et al. 
[19793 using ISEE-1 measurements on November 20, 
1979. However, this was a very turbulent period 
with extremely violent and irregular magnetopause 
motion. Heikkila [1982] has pointed out that 
these variations cause inductive effects in the 
data so that the interpretation of the 
measurements is not as simple as first believed. 

Ffest of the work on the energetic electron 
layer vriiich surrounds the magnetopause was done 
prior to the ISEE results [Bieber and Stone, 
1979; Meng, 19793. Scholer et al. [1982b] have 
examined high resolution energetic electron 
measurements for one ISEE pass through the 
subsolar magnetopause. They find a very 
filamentary structure with scale size of few 
electron gyro radii. Ion composition 
measurements have been performed in the bomdary 
layer, magnetopause and adjacent magnetosheath 
[Peterson et al., 1928]. All three regions 
contain both ionospheric and solar wind 
components. Plasma waves are intense near the 
magnetopause [Gurnett et al., 1979; Tsurutani et 
al., 1981; Anderson et al., 19823. In flux 
transfer events, the dominant plaana wave 
features are an intense low frequency continum, 
a dramatic increase in the ft*equency of 
occurrence of short wavelength spikes, 
quasi-periodic electron cyclotron harmonics 
correlated with "1 Hz magnetic field fluctuations 
and enhanced electron plasma oscillations 
[Anderson et al., 19823. 
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RECONNiCTION 


There has long been little doubt that 
reconnection as described by Dungey [1961] is an 
important process in the magnetosphere [cf. 
Russell, 1976]. However, until the launch of 
ISEE, the existing plasma data lack sufficient 
time resolution and 3-D coverage to monitor the 
variations in plaana behavior predicte in a 
reconnecting magnetopause. As expected ISEE^I 
and >2 soon encountered the expected signatures 
of merging or reconnection in the plasma data 
[Paschmann et al., 1979] as shown in Figure 9* 
However, the quasi-steady state merging signature 
was clearly found on only a few occasions in two 
seasons of observations [Sonnerup et al., 1981]. 
On the other hand, when conditions are right, the 
reconnection process can proceed continually for 
many hours [Gosling et al., 1982]. Further, some 
of the reconnection must be taking place in an 
unsteady manner in FTE's. On the other hand, 
Paparoastorakis et al. [1982] show that 
reconnection doesn’t always occur when you might 
expect it to. 

Eastman and Frank [1982] have questioned 
PaschDoann et al.'s [1979] identification of 
reconnection. Their critic ian in turn has been 
examined in detail by Scholar et al. [1982] and 
Daly and Fritz [1982], and successfully couitered 
by the later two authors. 

Despite the success of ISEE observations in 
demonstrating the reality of reconnecton, there 
still remains the problem of understanding how it 
operates and what controls it. This problem is 
being attacked not Just through observational 
programs but also through analytic theory 
[Coroniti, 1980; Grenly and Sonnerup, 1981; Quest 
and Coroniti, 1981a, b], computer modeling [Cheng, 
1979; Sato and Hasegawa, 1980; Birn and Hones, 
1981; Brecht et al. 1982; Matthaeus, 1982; Sato 
and Walker, 1982] and laboratory studies 
[Gekelman and Stenzel, 1981; Gtkelman et al., 
1982; Stenzel and Gekelman, 1981, Stenzel et al., 
1982a, b; Bam and Bratenahl, 1982]. The 
experiment of Bam and Bratenahl is directed more 
to dayside reconnection, tdiereas the experiment 
of Stenzel, Gekelman and co-workers is more 
appropriate to the magnetotail. The latter 
experiment has yielded a rich harvest of results 
Including the observation of double layers in the 
current sheet of the reconnecting jdasma. 

If reconnection occurs on the dayside, then 
reconnection must also occur on the night side of 
the magnetosphere also becuase dayside 
reconnection if left unchecked and unreplenished 
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would drive all the magnetic flux into the 
magnetotail. Observations in the tail assure us 
that indeed reconnection takes place there. Caan 
et al, [1979] I Hones and Schindler [1979]| Hones 
[1980], Nishid^a et al. [1981] and Hayakawa et al. 
[1982] observe the expected joint field and 
plasma behavior. Bieber and Stone [1980; 1982] 
Observe the expected streaming electrons and 
plaana behavior. Forbes et al. [1981a,b] observe 
the tailward retreat of the neutral line and the 
expected (daana sheet drift. There can be little 
doubt from the in situ tail data that 
reconnection is also occuring there. 

PUSNA SHEET AND NEUTRAL SHEET 

The magnetotail is extremely Important in 
magnetospheric physics as the site of energy 
storage for substonn related processes. Thus, 
some of the tail related papers will be discussed 
in the reviews of Burch [this issue, 1982] and 
Hughes [this issue, 1982]. However, some of 
these works merit the risk of repetition. 
Fairfield et al. [1981a] combined IKP-6 field and 
plasma data to study the accunulation and release 
of energy in the tail. Coroniti et al. [1980] 
combined IMP-7 plasma field, energetic particle 
nd plasma wave data to take a detailed look at 
plasma sheet behavior during substorms. Rapid 
flows and highly turbulent fields were observed. 
Erickson and Wolf [1980] questioned on 
theoretical ground whether steady-state 
reconnection is even possible in the tail. Hones 
et al. [1982] report observations of the 
three-dimensional plasmoid that is formed during 
reconnecton in the tail. 

The neutral sheet has long been known to be 
pulled north and south by the diurnal wobble of 
the earth's magnetic equator as if the neutral 
sheet were hinged to the magnetic equator at 
about 10 Re. Fairfield [1980] has used IMP tail 
observations to refine the model of how the 
neutral sheet responds to this diurnal torquing. 

Hardy et al. [1979a] have examined the 
plaana mantle and bomdary layer plasma as seen 
in the surface of the moon with the suprathermal 
ion detector experiment. They find that the 
appearance cf these ions is controlled by the 
interplanetary magnetic field. The probability 
of appearance is gt'eater for southward fields and 
in the northern lobe is greater on the dawn side 
of the tail when the solar magnetospheric 
Y-component of the IMF is positive. In the south 
a positive Y-component increases the probability 
in the dusk sector. Closer to the earth 
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Fairfield et a} . [1981] report that the plasoa 
sheet can become as thin as 1000 Ion during 
substorms. Ihe energetic ion composition of the 
plasma sheet has been investigated by Peterson et 
al. [1981]. Ihey find that the plasma sheet has 
a variable ionospheric component representing 
from 10% to more than 50% of the total nunber 
density and that there must be more than one 
process responsible for the energization of solar 
wind plaana to plaana sheet energies. 

A thin layer of earthward stresning 
energetic protons and alpha particles has been 
observed by several groups Just external to the 
plaana sheet [Hobius, 1980; Spjeldvik and Fritz, 
1981; Williams, I 98 I]. The layer a;^ars to be 
about two gyro radii thick and can be seen to 
reflect from the earth and stream tailward , 
Another important observation in this sane region 
of space is the detection of large electric field 
spikes [up to "SO mV o'*'] usually withir one 
minute of the plasma sheet boundary [Cattell et 
al., 1982]. These strong electric fields occur 
in regions of enhanced low frequency turbulence 
and in regions of field«aligned current. 
Individual electric field spikes are 
well-correlated with anall-scale gradients in 
particle fluxes and small-scale currents. 

Sharp et al. [1981] have examined low energy 
ion streams in the magnetotail boundary layer, 
lobes, and plasma sheet. They find that the 
boundary layer or plasma mantle consists of 
plasma of solar wind origin but that the streams 
in the lobe and plasma sheet have ionospheric 
composition. Finally, we note that Meng [1981] 
has presented a statistical survey of the 
energetic particle population of the magnetotail. 

OTHER PUSMA BOUNDARIES 

Other plasma bouKlaries that were the object 
of intensive investigation in the past received 
little attention over the last four years. Ihe 
polar cusp has been studied mainly with low 
altitude satellites and rockets. Atmospheric 
Explorer-D low energy electron and ion data have 
been exanined by Reiff et al. [I960] and Burch et 
al. [1980] to determine the effect of the 
inter^anetary magnetic field on the cusp. When 
the interplanetary magnetic field is southward 
the average energy of ions decreases towards the 
pole. When the interplanetary magnetic field is 
northward it decreases and then increases. The 
former signature is interpreted to be consistent 
with merging and the latter with diffusion. .Ion 
composition measurements from S3-3 show that at 
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tines the magnetosheath plasma enters 
adiabatically while on other occasions the ions 
appear to have passed through an electrostatic 
potential in reaching low altitudes [Shelley, 
1979]* These ion composition data, however, also 
show that ionospheric ions are continuously being 
accelerated to energies of the order of keV with 
the cusp region. Heng [1980] has examined the 
variation in polar cusp position during a 
geomagnetic storm. Curtis et al. [1979] have 
examined high altitude observations of plasma 
waves in the cusp during substorms. A region 
almost devoid of plasma has been disct^vered by 
Calvert [1981] at high invariant latitudes in the 
night magnetosphere on the basis of inferred 
density determined from natural plasma wave 
observations. This auroral plasma cavity, as it 
has been called, has densities below 1 cm~^ at 
distances of 2 Re and above and is believed to be 
a transient phenomenon associated with the 
generation of AJCR. 

Another topic that has recieved very little 
attention lately is the injection boundary of 
substorm particles and the nose events of 
Explorer This appears to be in part because 

it is now fairly well imderstood and in part 
because new and exciting plaars injection 

mechanisms have been discovered. In the 
steady-state electric and magnetic fields of the 
magneto s|;^ere plasma can drift from the plasma 
sheet and around outside of some demarcation 
boundary. Inside this bouidary which varies with 
pitch angle and energy, plasma circulates 

enclosed paths and does not intersect with the 
plaana sheet. During substorms the convection 
electric field in the magnetosphere increases and 
decreases on time scales comparable to the drift 
time of this plasma.. Kaye and Kivelson [1979]| 
Southwood and Kaye [1979] • Kivelson et al. 
[1983] and Ejiri et al. [1980] successfully use 
simple electric field models and time variations 
of these models to replicate the observed 
substorm-associated features of the low energy 
plasma. The more comprehensive Rice model is 
discussed by Walker [this issue, 1982]. 

In the area of understanding ion injection 
in the magnetosphere attention is now being 
focussed on field-aligned flows from the 

ionosphere. When these ion beans Mak at sane 
angle intermediate between 0^ and 90° pitch angle 
they are termed conics. Such distributions have 
been seen deep in the magnetosphere with the 
ISEE-1 plasma composition experiment in all three 
primary ionospheric species [Horwitz et al., 
1982]. Measurements with the ATS-6 plasma 
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spectrometer also reveal field-aligned thermal 
ions in the midnight region [Olsen, 19S?]« An 
even more exotic pitch angle distribution h:> been 
discovered with the ion analyzers on the P78-2 
satellite [a.k.a. SCATHA] and tei*med ion 
* zippers' [Fennell et al., 1981]. These ions are 
predominantly field-aligned at low energies and 
predominantly peaked perpendicular to the field 
at high energies with a very narrow transition in 
energy. They have been called zipper events 
because of their appearance oi. energy-time 
spectorgraus made using a detector that scans in 
pitch angle. The two components of the zipper 
are quite distinct in their magnetospberic drift 
paths. The low energy component drifts to the 
dayside via local morning and the high energy 
component by local evening. Ion composition 
measurements show that the low energy parallel 
component is mainly composed of oxygen and the 
high energy perpendicular component consists 
mainly of hydrogen. Lyons and Mbore [1981] 
attribute this difference to the effects of 
charge exchange. 


THE FUTURE 

While much has been learned about the outer 
magnetospberic boundaries from the ISEE-1 and -2 
missions, the analysis of data from these 
spacecraft has far to go. Thus, we should expect 
many new results still from this mission. The 
Dynamics Explorer mission has just been launche ' 
ana we should learn much auroral plaanas at mid 
and low altituies in the near future. However, 
we will not learn much about the distant polar 
casp until new spacecraft are launched . 
Hopefully, the Polar Plasma Laboratory of the 
OPEN mission which is now being planned will fill 
this void. 
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/ibstract. Corparisen of Poi.’er spectra cf 
nagiiotic fiela data fror IEEE-1 end -2 recorced 
sinultnneously on both sides cf tlv: riagnetopanse 
showed tl.at pDvier level inside thse nannetosphere 
varied v;ith pa.’er level outside in the macne- 
toshceth and suggested tJict tl:c- sai!C frecuencies 
v/erc erih^lnced or. the t\i?c sides of the boundary, 
Pov’cr levels v;cre tvo to three orders of rr.gni- 
tude la/er inside tlian outside the ragneto- 
sphcrc, indicating tlat v/ave energy was 
transritted inside froni the sheath. 


HTjT^CXXJCTIOIT 


A persistent and significant, although v;cak 
and disordered, correlation between solar wind 
properties and dayt.ine geonagnetic pulsaticxi 
activity in the Pc 3,4,5 ranee, periods T=10 to 
500sec, has been established by rnnny reports, of 
v;hich v;e cite a sarnpling iBol'sliakova and 
Troitskei’a, 19FC; Gul'elmi, 1974; VJett) and Orr, 
1976; Saito et al., 1979; Greenstedt et el., 
1979; Uolfe et al.. If 00]. Tliese correlations, 
together with tJ)e need to understand pulsaton 
phe.nonena in general, have also led to a ntnrber 
of models attempting to explain the observations 
ISouthv;ood, 1968; Greenstadt, 1972; Vinogradov 
and Parkhonov, 1974; Chen and Kasecavra, 1974; 
Kovner, 1976], Doth observation and theory have 
been concerned with the sources of the vraves, 
the modulation of the waves by the conplex media 
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through which they travel r and the control of 
wave pr(^'rtie& directly or indirectly by solar 
wind conditions. A nore conrreh€.nEive surr.ic.ry 
and reference list can be found in a paper by 
Greenstadt et al. [I960]. 


Pegardless of tlic specifics of either 
rodels or observations, ha^cver, there arc, in 
the most general terms, three possible sources 
of pulsations in the magnetosphere: Vfeives are 
produced inside the ragnetos^iere, waves enter 
the magneto^)here from the magnetosheatli , or 
waves er.icJiate from tliO bouridsry surface separate 
inc the magnctoshcath and the magnetosphere 
(i.c. from the magnetopause), ftone of these 
possible 50>jrces e>:cludes either of the others, 
but certainly if waves enter the magnetosphere 
free, outside, they must cross the boundary some- 
wiiere at some time, and it should be possible to 
c:;tuhlish by obsc-rvaticT! a physical sirilarity 
betwe-en waves detected siriiltantously on botl. 
sides of the magnetopause. This is the approach 
we pursue in the proser;t study; the preliminary 
results reported here support the transfer of 
v.-ave energy im-’ard frer. the rr-gnetosheath to the 
nagneto^here. 

Ue define as a "straddle" a situction in 
which one spacecraft is on one side of the meg- 
netopause and a second is or» tlie opposite side. 
71 le ISFC-l and -2 setellitos providc<] neny such 
straddle crossings, but ve have concentrate on 
th:e magnetometer data of 1970, wt»en the space- 
craft \:ere separated by hundreds to tliousands of 
km for several montlis, giving straddles long 
enough for unequivocal analyses of waves with 
perils up to several minutes. Our approach is 
to compute concurrent pa-’er spc-ctra on both 
sides of clearly defined magnetopause crossings 
and to seek siralarities or differences in the 
^jectra from the two spacecraft, \ie require rea- 
sonably .continuous data frem both satellites 
over minimal intervals of 15 minutes and prefer- 
ably ruch longer. V)e found four potentially 
suitable stra^les of which we have developed 
two for this preliminary report, supplemented by 
two singlc-^cecraft crossings, as explained 
later, Uc present only spectra of the total 
field magnitude at this time, derived from vec- 
tor saiiples every two seconds, 

DATA 
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The tvt> cases we describe occurred on 6 
October and 27 Hoverticr 1978. The 8 October 
ctossino occurred in early afternoon, the 27 
Moveirber crossing in lete morning. The da^-time 
surface field for our cases vras measured b^' the 
AFGL ground stations in a sector ipanninc the 
Lhited States near 55 degrees iregnctic latitude. 
CrourK'i station data processed for this report 
v;ere recorded bj' the Nevport station at the 
western edge of the sector, vihich fell before 
and after noon during the two intervals of 
interest, so that the morning-afternoon local 
times on the ground were reversed from those of 
our crossings at the nacnctopausc. Vte used 
ground station data as a rough guide to the 
presence of traditional pulsations on the 
ground, since it is v;aves related to the origin 
of such ixilsotior^ that wish to be studying. 
T^ie two straddle intervals of this report to^ 
place during ITT conditions favorable to the 
presence of enhanced v/ave activity in the subso- 
lar magnet osheeth convccted frer. quasi-perollel 
bov; shock structure. 

Pcn;cr spectra presented in this study were 
calculated with the fast Fourier transform. 
Figure 1 places one spcctrur.., dotted curv^e, ot 
the sar.e scales use<l by Fairfield fl97G] to 
display the characteristics of magnetosheath 
rrrnetic noise from several spacecraft. The 
dotted curve is a po\;er spectrum for tlie total 
ambient field recorded by ISEE-1 irmcdiately 
outside the magnetopause on 8 October. All 
spectra are for the total field magnitude; 
details of the earlier spectra can be found in 
Fairfield's reviev; [op.cit.l. 

Hie figure illustrates four principal 
points: magnetosheath spectra topically show 
either an enhancertent or a slope change, or 
bC'th, at or below the local proton cyrcfreguency 
(highlighted b\' the dashed lines) ; second, mac- 
netosheath spectra arc highly variable, both in 
absolute pov;er level and in the freouencies that 
right be cnhcinced in any particular sar.ple; 
third, spectra taken from ISTH data >’ears later 
are reasonalDly representative of the same v/ave 
beliavior tliat prevailed during the earlier meas- 
urements. Fourth, and most importantly for this 
report, spectra obtained from ISEE close to tlie 
magnetopause do not appear to define a special 
region in any v/ay unrepresentative of the magne- 
tosheath at other locations. 
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Piaure 2a is an escesnple of our first strad- 
dle case. The upper panels of the figure display 
plots of field ragnitudc frcr.! ISEF-1 and ISTE-2, 
for the magnetopause crossing of 8 October 1978. 
ISCE-2r loiTer field plot, entered the racneto- 
snhere first at 1805t50 and finally at 1813; 
lSEE-1, top, encountered the magnetosphere 
later, initially at 1833 and entered finally at 
3635:40. Tlius, there were 38 minutes during 
which data were acquired sirtultancoucly frcci one 
satellite outside and one inside the magneto- 
pause. 


Spectrur. h chovr. the wa\»e pov.'er in the 
total field in the racnetoeheatli just outside 
the magnetopause, at ISEE-1. The next spectrum 
bclov;,' B, shows the wave pa;er in the mgneto- 
sp^re just inside the macnetopause, at ISCE-2, 
for the same time interval as that of the first 
spectrur.. ^le pov:er was cppreciaLly la-tr and 
the decrease in pov.*er with frequency clearly 
much steeper inside th^n outside the magneto- 
pause, twinning with about one third the out- 
side power at tlie lowest frequencies. At 0.1 Hi., 
there v;crc throe orders of magnitude difference 
between the two spectra. 


Spectrum G, at bottom, represents the power 
on t.he ground at the TSCL station at llex.-port, 
’.’ashincton, for the sane intervals as in tJie 
dcpict^'i satellite samples. The ground station 
v;as a fev; hours v.»cst of the satellites, about 
local noon. !Ihe shaded vertical strips in the 
spectral panels drav? attei'.tion to tlic enhance- 
ments in pov;er at the satellites and on the 
ground that appear to bind vratve activity 
together in the inner nacnetosheath , the outer- 
most magnetosptjere, and at the earth's surface. 
All the spectra shw some concentration of power 
between .02 and .07 Kz in the form of a plateau 
or peaks in the respective curves. 

Figure 2a shov.*s a progressive decline in 
po\rcr from the nagnetosheath to the earth's sur- 
face. Each spectral curve is contained in, i.e. 
accounts for a fraction of the pcx/er of, the 
next spectrum above it. The magnetospheric spec- 
tra are well below that of the nagnetosheath and 
are closer to one another than to the latter, 
the discrepancy being greatest at the highest 
frequencies. 

Figure 2b superposes spectra from a second 
straddle case on 27 Moverrber 1978, v4)en ISEE-1 
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was just outside the nagnetocause. ISEE-2 «os 
deeper inside the nagnetos^^rer having entered 
at 2000 tl)£n it had been on 8 Octw and thie 
Newport station was below and east of the satel- 
lite meridian (in the early afternoon sector). 
In this instance, the pm;er in the sheath (A) 
displayed enhancement and a plateau between .011 
and .05 Ilz, as did also the pex^r in the ranne- 
tosphere (B) , v.t»ile the cot responding pov.»er’ on 
the ground (G) v/as relctivily featureless, but 
essentially at thie same le\*cl as at IFFT-2. 
l?het)ier the apparent lack of frequency' enhance- 
ment on the grouTid at tliis time was because of a 
dele>*ed effect not yet visible, an unfavorable 
position in the afternoon sector, or a poor 
choice of representation of the surface record 
is still to be dotemiined. Tlie small grap^ at 
the bottom shows the power distribution in Ej* at 
I■e^port for the local noon and afternoon inter- 
val including t)»e 22-rinute sepront cf the upper 
gra0»; clearly, tliere was some activity in tJic 
surface field within the longer interval and 
v;itl:in the enhanced portion of the spectrum at 
the satellites. Ihe attenuation of wave pa;cr 
across the nernctoprusc is ob\»ious here, as in 
the previous case, but ws also see that the fre- 
queno' rar.ge of erhex.cerent v;as shifted to some— 
v;hut IcxTcr frc<;ucncics, in all locations, th.an 
in thi6 case of C October, ao indicated bj' the 
shading in Figure 2. 

The megnetosheatlj spectra of B Oct. and 27 
l%)v. ere superposed in Fig. 3, showing that the 
power of tiic- 27 Itov. spectrum peaked at lower 
frequency and dropped more rapidely thari the 
power cn 8 Oct. This difference corresponds to 
the slightly different regions of the frequency’ 
scale that seemed to be shw; enhanced pov.’er in 
the magnor^Jhere, as shaded in Figure 2 . Unfor- 
tunately, the absolute power levels in the irog- 
netosheath in the two cases exairiinod above did 
not differ appreciably from each other, consid- 
ering the wide range of power exhibited in the 
curves of Figure 1. In order to study whether 
poi;er inside the magnetorphere is related gen- 
erally to pov;cr in t)>o magnctoslKath, it was 
necessary to exardne nonstraddle cases. Oi the 
premise that tlie sheath spectrum remained sub- 
stantially unchianged from one interval to the 
next over an hour's time, we selected cases with 
significantly differert paver levels and com- 
pered spectra before and after magnetopause 
crossing, rather than 8ir.vltancous spectra on 
opposite sides of the boundary. Corresponding 
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magnetoeheath and magnetosphere ^sectra are 
superposed in Figure 4 for three days, 6 
October, and 10 and 17 Septcnt>er 1976. The 6 
October curves are already familiar. The new 
ones shok* that progressively lo%;er pov?er in the 
sheath corresponded to progressively la;cr TX>wer 
in the magnetosphere, suggesting, with these few 
cases and the necessary assunption of stationary 
ity on the lOtli and 17th, that the pov«rs inside 
and outside the magnetopause were directly 
related. 




The data presented above may be sumarized 
as follows, with the understanding that we refer 
essentially to the frequency range 0.01 < f < 
0.1 Hz (periods 10 < T < 100 sec.): 

The power within the magnetocheatl) was 
10 to 1000 times the pov/er in the mag- 
netosphere; 

The po'.;cr within tlic r^sgretospliere varied 
less then a factor of 10 fror.i the magneto- 
pause to the surface; 

Tlie pov;er level inside the macnetosr^jere 
correlated overrl3 with pcn;er level outside 
the magnetopause; 


•Rje frequency of power enhancement in the 
nacnetosphere appeared to shift with the 
frequency' of pov;cr enhancement in the mag- 
netoshcath. 

Tl^ poiTcr outside the mnenetopause appered 
to be representative of pa«r in the megne- 
tosheoth generally. 

In addition to the- foregoing, we have found evi- 
dence tliat the variable presence of e frequency range 
of cnhfjtced pov^jr in the mngnetosheath was corrclat^i 
with the variable presence of an EIF orientation 
favorable to tlic occurrence of tjuasi-parallcl struc- 
ture in the bow shod; around the subsolar point. This 
result will be treated in a separate report. 

DISCUSSIOJ 


Similarity between ^>ectra in the magne- 
tOGheath and magnetos[4)ere may be explained. 
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excluding coincidence, in the three ymyt defined 
in the IinrCDUCTICl). observations listed in 
tlie SUtl'Anv sucger.t tliat internal nagnctcsplteric 
origin nay be excluded from consideration 
because of bie higher power observed outside the 
boundary, lioreover, a principal criterion by 
which wc selected cases was the appearance of a 
clear nagnotoneusc allowing us easily to deter- 
rine that one ^cecraft was inside, the other 
outside the boundary. This would eliirdnatc wave 
propagation outward along field lines locally 
interconnecting the solar wind to the ntagneto* 
splje-re. In fact, wc have established that in 
one case the boundary' was well aeproximated by a 
tangential discontinuity. 

Surface waves are an unlikely explaration 
because of the overv/J-ielring pa/er in the nacne- 
tosh«eath comared to the magnetosphere. The 
nc?\;e£t calculations cf surface v.’sve effects [Pu 
and r.ivelson, 19C3] retire that tlie magnetic 
v.’avc pa.’er inside- the boundary' exceed that out- 
side, opposite our results. Also, the freouen- 
cies e::pcctec for surface vraves tend to be la%»er 
than those with v.hich \k have been dealing. 

I’e conclude therefore that our preliminary 
results are consistent \;ith eiiternal vave ori- 
gin, specifically with the transfer of a small 
fraction of rr«icnetosheath wave pcn-.*er, possibly 
derived from quasi-pc.rallel shock structure, 
into the magnetosphere to appear as waves in the 
Pc 3-4 range. 

Tlie asserted cormonality of frequency 
enhancesnent across the boundary, as illustrated 
here in Figure 2, is subtle at best. This is 
hardly surprising, since we are dealing with c 
global phenomenon notoriously elusive to sharply 
^fined correlations, vMch v;c chance to sarple 
as a fev; straddles ct a few points in qiace. V?e 
chose the two examples here as the purest strad- 
dle cases. Other spectra, with more persuasive 
enhiancement profiles, were obtained in data con- 
texts requiring r.torc e>:position then could bo 
included in Uiis letter and will be the subject 
of a separate report. The component-bj'- 
component details of the transfer process, the 
global picture describing v^re the most effec- 
tive transfer takes place, and the pathwa^^ 
whereby' broadband energy in the magnetosheath is 
recorded as monochromatic pulsatio^ in the mag- 
nctos^iere rer.»in to be determined. 
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Figure 1. /in ISSE-1 racnetorhesth spectrum, 
obtairifcd near tiiC ^^.gnctorzuse (dots) , super- 
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obtained near the magnetopause (dots) 
posed on a selection af spectra from 
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Figure 2. Magnetic field magnitude records and 
superposed Q>ectra for two jcases of ISCE~1» 
ISEE-2 straddles of the magnetopause x (a) 6 
Oct.r (b) 27 Nov.r 1978; k, B, and G signify, 

qpectra for the indicated intervals at ISEE-l, 
25CE>2, and the APGL Newport ground station. The 
insert at bottom right covers a longer interval 
at Newport f as noted. 
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The spacecraft ISEE-3 will be sent to encounter G1acob1n1>Zinner 
In September 1985 during the period of Intense Halley observations. This 
mission provides a unique opportunity for comparative study of Important 
dynamical processes at both short and 1ong-pev*1od comets. The ISEE-3 
payload Includes a vector Helium magnetometer and a high-sensitivlty plasma 
wave search coll that give complete Information on the magnetic field 
profile and the B-fleld fluctuation spectrum up to 1 kHz. Long electric 
dipoles (90 meters, t1p-to*t1p) are used by the plasma wave Instrument and 
the radio mapping receiver to detect wave phenomena over the spectral 
range 17 Hz to 2 MHz, and the electron section of the ISEE-3 plasma probe 
measures the plasma density, the flow speed, and other Important character- 
istics of the electron distribution function. These Instruments can 
readily Identify the position and strength of the bow shock, provide data 
on the nature of lonosheath/coma flow and turbulence, and yield direct 
Information on the plasma phenomena that develop at the contact surface 
and In the foreshock region. As ISEE-3 flies across the comet tall. 
Identification of the upper hybrid resonance emissions or trapped continuum 
radiation by the two wave Instruments should provide unambiguous and 
absolute data on the local electron density profile, similar to the science 
return from the Voyager 1 wave Instruments In the wake and tall of Titan. 
Other ISEE-3 Instruments that measure characteristics of energetic particles 
and Ion composition will also yield very significant Information on the 
solar wind-comet Interaction. 

ISEE-3 (the third In the series of International Sun-Earth 
Explorers) was launched on August 12, 1978, and Figure 1 contains a 
drawing of the soacecraft with all booms and antenna systems deployed. 

Figure 1 also has a drawing of Comet Glacoblnl-ZInner. This sketch 
Is based on photographs taken by E. Roemer during the 1959 apparition 
(see Figure 6) and It shows that this active short-period comet has 
many very Interesting characteristics Including an extensive visible 
coma (outer diameter z 50,000 km) and a long, narrow Ion tall (length 
of visible portion * 350,000 km at 1 AU). 

The ISEE-3 encounter with Glacoblnl-ZInner will represent the 
third extended mission phase for the spacecraft. During Its primary 
mission (August 1978 through August 1981) the spacecraft was stationed 
In an elliptical orbit around the Li llbratlon point (approximately 
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Figure 1. ISEE-3 and Comet. Giacob ini-Zinner 


250 Re upstream from earth), and from this upstream position it monitored 
bolar wind onenomena and processes and provided early warning of inter- 
planetary storms. During the first phase of the extended mission it 
remained in the upstream halo orbit, but in June 1982 the thrusters were 
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fired, sending ISEE-3 back to the earth's geomagnetic tall. The first 
tall crossing beyond the lunar orbit occurred during October 1982, and 
ISEE-3 will continue to explore the distant geomagnetic tall out to L 2 
(250 Re downstream) until Deceiiter 1983, when a close lunar flyby will 
launch the spacecraft Into Its mission phase toward Comet Glacoblnl-ZInner. 
The trajectory for this phase Is shown In Figure 2. Here the coordinate 
system Is earth -centered and sun-oriented, and It can be seen that when 
ISEE-3 encounters the comet on September 11, 1985 (six days after Glacoblnl 
Zinner perihelion), the viewing from earth will be excellent. Figure 2 
also shows that ISEE-3 should be able to provide significant Information 
on solar wind conditions upstream from Halley's comet after the Glacoblnl- 
ZInner encounter; the Halley-related observations will represent a fourth 
extended mission phase. 
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Figure 2. ISEE-3 Trajectory Relative to Fixed Sun-Earth Line 
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The Glacoblnl-ZInner encounter will occur as the comet crosses the 
ecliptic plane (descending node) at a heliocentric distance of 1.03 AU 
and a geocentric distance of 0.46 AU. In the frame of reference of the 
comet, the spacecraft will be traveling at a speed of 20.7 km/sec, primarily 
from south to north. It Is planned to target ISEE-3 to pass downstream 
of the nucleus, and Figure 3 shows two views of the planned encounter. We 
expect that ISEE-3 will have Inbound and outbound crossings of the bow 
shock, the contact surface, and the boundaries of the outer coma or Inner 
plasma sheet. In addition, ISEE-3 will traverse the foreshock, the sub- 
sonic sheath regions, and the tall lobes. ISEE-3 will be targeted to fly 
by at a downstream distance which Is deemed to be safe In terms of the 
dust hazard. A distance of 10,000 km has been tentatively proposed, but 
the dust hazard has not yet been carefully evaluated. Nevertheless, It 
Is clear that for Glacoblnl-ZInner, the dust problem Is considerably less 
than the corresponding problem at Halley. For Glacoblnl-ZInner, the flyby 
speed Is 3.5 to 4 times lower, the rate of production of dust Is perhaps 
two orders of magnitude less, and observations of the Glacoblnid meteor 
stream suggest that dust particles from Glacoblnl-ZInner may have a mean 
mass density two orders of magnitude lower than the average for the Halley 
dust particles. 
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Figure 3. Comet Intercept 
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In Table 1 the science payload of ISEE-3 Is listed and the principal 
InvestiMtors are Identified. For all but the four cosmic ray Instruments 
(bottom), the functional measurements are described, and the anticipated 
encounter sample times and spatial resolution scales (for 1024 bit/sec 
telemetry rate and 21 km/sec flyby speed) are tabulated (the Ogllvie 
Instrument will be configured to operate In Its most rapid scan mode). 

This table Indicates that ISEE-3 has an excellent science capability for 
study of the solar wind-comet Interaction, and we proceed to Illustrate 
this capability by showing In Figure 4 some simultaneous ISEE-3 measurements 
obtained during a 30-hour Interval on November 11-12. 1978. when an Inter- 
planetary shock swept past the spacecraft. These Interplanetary measure- 
ments resemble those that will be obtained In the foreshock, near the shock, 
and downstream. 


TABLE 1; ISEE-3 SCIENCE PAYLOAD 


Investigation 

Principal 

Investigator 

1024 BPS Sampling 
Period (sec) 

Spatial 

Resolution (km) 

PLASMA ELECTRONS 
5 EV to 1.0 keV 

S. J. Bam* 

24 SEC 

500 KM 

MAGNETIC FIELD 

1 M GAMHA UP TO t 1.4 GAUSS; 0-3 HZ 

E. J. Smith 

1/3 SEC 

7 KM 

PLASMA WAVES 

E-FieLD; 16 CHANNELS; 18 HZ TO 100 kHZ 
B-FIELO: 11 CHANNELS; 0.3 HZ TO 1 kHZ 

F. L. Scarf 

1 SEC 

21 KM 

PAOIO WAVES 

E'FielO: 24 CHANNELS; 30 kHz to 2 MHZ 

J. L. Steinberg 

112 SEC 

2.400 KM 

PLASMA IONS 

Velocity; 20-200 km/sec ; M/Q: 4-50 

K. w. Ogilvie 

1.200 SEC 

25.000 KM 

ENERGETIC PROTONS 
30 KEV TO 1.4 riEV 

R. J. Hynds 

32 SEC 

660 KM 

ENERGETIC ION COMPOSITION 
30 KEV/O < E/0 < 600 KEV/0 

D. Hovestaot 

120 SEC 

2.500 KM 

X RAYS AND ELECTRONS 
Electrons: > 200 keV 
X Rays; $-1250 keV 

K. A. Anderson 

1 SEC 

21 KM 

HIGH ENERGY COSMIC RAYS 

E. C. Stone 



HIGH ENERGY COSMIC RAYS 

H. H, Heckman 



MEDIUM ENERGY COSMIC RAYS 

T. VON Rosemvinge 



COSMIC RAY ELECTRONS 

P. Meyer 
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Figure 4. ISEE-3 Measurements During Passage of an Interplanetary Shock 


The magnetic field profile in the bottom panel of Figure 4 shows 
that the interplanetary shock arrived at ISEE-3 at 0028 on November 12, 
and that the extensive downstream region of disturbance was characterized 
by the presence of intense low-frcnuency magnetic turbulence. The velocity 
and density panels from the plasma electron instrument also indicate 
passage of a strong shock with post-shock flow disturbances, and the 
18 Hz panel from the plasma wave instrument shows that the region had 
large-amplitude magnetic noise that we identify as whistler mode waves. 
Other significant post-shock characteristics are the elevated 60 kHz 
noise levels from the radio mapping investigation and the very large 
enhancements in the fluxes of energetic protons and helium ions. In 
Figure 4 there is also clear evidence of a foreshock region with high 
amplitudes for the ion acoustic waves (3 kHz E-field) and a pre-shock 
build-up for the 30 keV protons suggesting local acceleration. The ion 
mass spectrometer profiles on the right side of Figure 4 illustrate 
the large variations in flow speed and composition over the two-day 
interval; at this time the instrument was dearly able to resolve changes 
in the populations ot ^He"'’'’', ‘*He++, O’’*’, and Fe ions. 
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The relevance of these interplanetary p1asn?a physics measurements 
for the Glacoblnl-ZInner encounter can be understood by considering the 
right side of Figure 5. Here ue depict the anticipated magnetic field 
and plasma density profiles by sketching possible variation-* of the 
electron plasma frequency, f„ [*(9i40 kHz, with N In cm"*], the electron 
cyclotron frequency, fr L*(2o B) Hz, with B In gamma], and related Ion 
and hybrid characteristic frequencies. Many models of the solar wind- 
comet Interaction suggest a shock-to-plasma sheet-to-shock distance of 
2-4 million kilometers for the downstream trajectory shown In Figure 2, 
and this would mean that the first half passage from the "shock region" 
to the center of the plasma sheet might take from 13-26 hours at a speed 
of 21.7 km/sec. Thus, the total traversal across the tall Is expected 
to occupy at least a full day. Moreover, the extensive foreshock contains 
comet-related energetic Ions, electrons, and enhanced turbulence levels 
for mhd and plasma waves, so that the total time for passage through 
the comet Interaction region may well Involve several days (The Space 
Science Board of the US National Academy of Sciences suggested that 
local measurements of the comet Interaction be made over a distance of 
10^ km. Implying coverage for five days at the ISEE-3 encounter speed. 
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Figure 5. ISEE-3 Wave Coverage (Left Side) and Anticipated Variation of 
Characteristic Wave Frequencies (Right Side) 


The labels on the right side of Figure 5 designate wave phenomena 
that must be studied during the comet flyby because (a) these measurements 
yield Information on very significant local wave-particle Interactions 
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(e.g., upstream waves; bow shock and lonosheath turbulence; current-driven 
plasina oscillations near the contact surface and the plasma sheet 
boundaries), or (b) because these measurements provide unique diagnostic 
Information (e.g.. detection of trapped continuum radiation In the tails 
of Earth and Jupiter yields absolute plasma density profiles; detection 
of upper hybrid emissions behind Titan yields similar absolute density 
Information). In addition, the very-low-frequency Ion and hybrid mode 
waves may be related to processes that develop continuously near the 
tall boundary region (plasma mantle) or Intermittently during substorm-1 Ike 
tall disruptions. Our experience In studying plasma physics phenomena 
behind Earth, Jupiter, and Titan suggests that comprehensive wave coverage 
Is needed, and the panels on the right side of Figure 5 show that the 
ISEE-3 payload will Indeed provide such coverage over the range from DC 
to 2 MHz. 

The expected science return from the ISEE-3 encounter with Comet 
Glacoblnl-ZInner Is sumnarlzed In Table 2. As noted above, ISEE-3 has no 
remote sensing or Imaging capability and no way to analyze physical or 
chemical characteristics of dust (although the plasma wave Instrument may 
be able to provide dust Impact data Just as the Voyager 2 Instrument 
provided similar Information during the crossing of Saturn's 6-r1ng). 
However, as shown In Table 2, the ISEE-3 Instruments will yield comprehensi 
Information on the dynamical plasma processes that control the solar wind 
Interaction with an active short-period comet. 


TABLE 2: EXPECTED SCIENCE AETURN - ISEE-3 TO GIACOBINI-ZINNER 

• Physics op the Cohetakv Bow Smock 

EvALUATI shock location, STRINGTH, STRUCTUNI, AND RELATED PARTICLE ACCELERATION; 

Determine characteristics or upstream ionization and mass loaoins; 

Determine extent or the poreshock. 

• Physics of the Iowosheath and Contact SunrAcs 

Determine the nature or the interaction ietheen the interplanetary magnetic pield 
AND the cometary PLASMA; 

Determine location, thickness, and microstructure op the inner iounoary (contact surpace); 
Evaluate sheath ionization phenomena and turiulence characteristics; 

Evaluate electron distribution punctions; search por electron peaks that might relate 
TO ionization OF PARENT MOLECULES. 

• Physics op the Comet Tail 

Determine configuration op the magnetic tail and content op the plasma sheet pilaments; 
Evaluate currents, have-particle interactions and acceleration processes in the tail; 
Evaluate ionic composition, energetic particle spectra, and electron heat plux; 

Search for plasma mantle regions, cross-tail E-pield acceleration processes, and 
substorm-like phenomena. 
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The encounter of ISEE-3 with Glacoblnl-ZInner creates opportunities 
for coordinated measurements that will enhance the overall scientific yield 
from the many 1985-86 comet Investigations now planned. Conditions for 
observing Glacoblnl-ZInner from Earth will be excellent during 1985. The 
bottom part of Figure 6 has one of the 1959 pictures of the comet* and 
the top of Figure 6 shows the most relevant parameters related to Earth- 
based viewing of the comet along with the corresponding Information for 
Halley. For an observer located at a latitude of 35N, Glacoblnl-ZInner 
can be viewed for four or more hours throughout 1985 and. In particular, 
at the tiine of the encounter. The predicted brightness (apparent total 
magnitude) will exceed that of Halley until October-November and will 
reach 8th magnitude at encounter. 

The phasing of Glacoblnl-ZInner relative to Halley Is such that 
viewing of the two comets tends to be mutually complementary. Glacoblnl- 
ZInner will be available for viewing during the dark hours for several 
months prior to Halley and for a month or so around the time of Halley 
perihelion passage which occurs In daylight. 

There will also be post-encounter opportunities for ISEE-3 to make 
particle and field observations In the solar wind upstream of Halley while 
Earth-based or space observations of Halley are being made. Figure 2 (above) 
shows the location of Halley during two Intervals of closest approach to 
the Halley-sun line. The first of these "radial" lineups would occur on 
31 October 1985 before Halley reaches perihelion, when the view of Halley 
from Earth will be excellent. At that time, the radial separation will be 
0.9 AU and fast streams, sector boundaries and other solar wind phenomena 
would be expected to reach Halley about 4 days later. The second lineup 
would occur 31 March 1986 after the Halley perihelion passage. This line- 
up will have a separation distance of 0.2 AU, and earlier In March 1986, 
as the other spacecraft encounter Halley's comet, ISEE-3 should be able to 
provide excellent upstream Information on corotating solar wind streams 
and transients. 
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•available dark hours - IS defined as the time interval during which 

THE SUN IS BELOW THE LOCAL HORIZON dV AT lEAST 18 DEGREES AND THE 
COMET IS SIMULTANEOUSLY ABOVE THE LOCAL HORIZON 



Mgure 6. Top: 1985-86 Observing Conditions for Comets Halley and 

Giacobini -Zinner . Bottom: 1969 Photograph of Comet 

Giacobini-Zinner. 
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ABSTRACT 

A comprehsnsive study of the November 11/12, 1978 shock event 
based on energetic particle,* solar wind, magnetic field and wave 
data from the ISEE-3, -1 and -2 spacecraft been undertaken 
both from the energetic particle and the collisionless shock 
point of view. The energy density of 10-50 keV protons acceler- 
ated by the shock Is found to be equivalent to the upstream mag- 
netic field energy density. The observations are in quantitative 
agreement with Lee's (1983) self-consistent theory for the excit- 
ation of hydromagnetic waves and the acceleration of ions upstream 
of Interplanetary shocks. 

1. Introduction 

We present selected results of a comprehensive study of the quasi - 
parallel interplanetary shock event of November 11/12, 1978 from both the 
collisionless shock and energetic particle points of view using magnetic 
field, plasma, MHD wave c energetic ion measurements on ISEE-3 ('v-200 Rc 
upstream from the Earth i 5-jw shock) and ISEE-1/2 (upstream, but close to 
the bow shock) (Kennel et al., 1983). Specifically we compare the obser- 
vations with quantitative predictions of current theories of particle 
acceleration by quasi-parallel shocks. Several papers have addressed the 
November 11/12 shock event (Reinhard et al., 1981; Kennel et al., 1982; 
Scholer et al., 1983; Tsurutani et al., 1983; Lee, 1983). 

2. Event Overview 

Figure 1 presents selected ISEE-3 energetic ion, plasma wave, solar 
wind and Interplanetary magnetic field (IMF) data for the period 1800 UT 
on Nov. 11, 1978 to 2400 UT on Nov. 12. The shock encounters lSEE-3 at 
0028 UT on Nov. 12. In the upstream period the 35-56 keV proton and the 
3 kHz ion acoustic wave Intensities Increased and the solar wind density 
decreased, all three with superimposed fine structure, until the shock. 
The IMF magnitude and solar wind speed remained virtually constant at 7.5 
nT and 375 km s*^, respectively. At the shock the Ion acoustic wave 
amplitude maximised and then suddenly dropped. The 35-56 keV proton flux 
rose by a factor of 10 In the last half hour before the shock passage to 
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Fig. 1. (left) 30-hr period of selected ISEE-3 data around the 
interplanetary shock encountered at 0028 UT on Nov. 12, 1978. 

N is the solar wind pseudo-density taken from the data pool. 

Fig. 2. (right) Time evolution for the period 2100 UT to 0100 UT 
of ISEE-3 energetic proton fluxes in 8 particle channels between 
35 and 1600 keV ((1) 35-56 keV; (2) 56-91 keV; (3) 91-147 keV; 

(4) 147-237 keV, etc.); ratio of differential fluxes in success- 
ive energy channels; first-order parallel anisotropy amplitudes 
of 35-56 keV proton distribution in solar wind frame; spectral 
density of 3.16 kHz ion acoustic waves; variance In the 0.01 - 
3 Hz IMF Y-component. 


a first maximum at the shock, reached a second peak about 0240 UT, comm- 
enced a slow decline and dropped suddenly near 1715 UT coincident with 
the arrival of the flare ejecta at ISEE-3. 

Figure 2 compares intensities, spectral indices and the anisotropies 
of the 35-1600 keV protons with ion acoustic and <3 Hz MHD wave ampli- 
tudes in the ISEE-3 upstream region. This entire region was free from 
effects associated with the bow shock. Detailed anal^«is of the solar 
wind electron observations, which exhibited strong bidirectional stream- 
ing, provides strong evidence that lSEE-3 entered a closed magnetic loop 
("bubble") at 2230 UT that intersected the approaching interplanetary 
shock at both ends. Inside this upstream bubble the ion acoustic and 
HKD wave amplitudes are intense and nearly constant till ^<0000 UT. The 
low energy (< 200 keV) proton channels respond significantly to the 







ORIGINAL PAG£ IS 
OF POOR QUALITY 


arrival of the bubble. On entering It the first-order proton anisotropy 
goes essentially to zero. Indicating an Isotropic distribution. In the 
last half hour prior to the shock encounter the Ion acoustic wave anpil- 
tude Increases, strongly correlated with the steep rise of the low- 
energy proton Intensity and the start of significant field-aligned 
streaming In the solar wind frame from the direction of the approaching 
shock. The <3 Hz MHD wave amplitude remains nearly constant. 




3. Shock Structure and Uave Activity 


Detailed analysis of the IMF and solar wind data from ISEE-3 and -1 
using four different technioues (Kennel et a1., 1983) resulted In the 
shock normal angle « 40 , a fast Mach number Mp » 2.6, a Jump of field 
strength of 2.2. Ue note that these parameters are for the local plasma 
subshock and do not take Into account the effects of energetic Ions. 


Ue suggest that energetic Ions must play a significant role In the 
shock structure. The measured 35-56 keV proton flux of the shock Is 
3 X 10** p/cn^ s sr. The corresponding number and energy densities are 
p cm*3 and 0.7 x 10“^® ergs cm*^, respectively. This energetic 
proton energy density, probably even an under-estimate. Is about 1/3 of 
the energy density of the upstream magnetic field, 2 x 10”^® erg cm”^. 
Extrapolating the observed proton spectrum down to 10 keV results In an 
increase of the particle energy density by a factor 4. 


Upstream waves are continuously present to unusually large distances 
from the shock. Emissions at 1900 UT on Nov. 11 imply that the scale of 
the upstream wave region corresponds to *v0.04 AU (Tsurutani et al., 1983). 
This region of enchanced upstream wave activity coincides with the ext- 
ended region of the energetic protons (Fig. 1). Figure 3 shows 8 min- 


utes of high-resolution magnetic field 
data surrounding the ISEE-3 (top) and 
ISEE-1 (bottom) shock encounters. 

Large amplitude (6 nT p-p) slightly 
compress ional waves with 'wT-s periods 
are apparent just ahead of the shock. 

The region between the shock and a 
strong IMF rotation contains extremely 
large ('^<20 nT p-p) transverse magnetic 
field fluctuations, with wavelengths 
(^•3 to 20 s periods) comparable with 
the thermal ion Larmor radius in the 
shock-heated plasma and therefore 
consistent with the resonant firehose 
Instability. The upstream and down- 
stream waves near the shock propagate 
essentially parallel to the local 
magnetic field. 






4. Comparison of Observations with 
Theory 

Lee (1983) recently presented a 
self-consistent theory for the excit- 
ation of hydromagnetlc waves and the 


Fig. 3. IMF Y-component and mag- 
nitude around the ISEE-3 (0028:20 
UT) and ISEE-1 (0058:30 UT) shock 
encounters. Both magnitude pro- 
files are essentially Identical. 
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acceleration of Ions upstream of Interplanetary shocks by first-order 
Fermi processes. We compare several facets of our observations with the 
predictions made by the equilibrium theory. 

a) The measured proton distribution transformed Into the solar wind 
frame shows field-aligned streaming away from the shock with an amplitude 
of *^<0.35 for the 35-56 keV protons In the ^100 R£ wide upstream region 
adjacent to the shock. In the near-shock downstream region we observe 
isotropy. If we transform our observations Into a frame moving with the 
shock velocity, we find Isotropy In the near-shock upstream region. 

These results agree with the prediction of Lee (1983) and the observa- 
tions of Scholar et a1. (1983). Farther upstream (2230-0000 UT) we also 
observe Isotropy in the solar wind frame. 

b) The energy spectrum can be fitted neither up nor downstream by a 
single exponent power law over the full 35-1600 keV range. For E t 250 
keV the differential proton spectrum has a spectral Index of y * 1.1 at 
the shocK and hardens with increasing upstream distance from the shock, 
in agreement ^fith theory. Tarther upstream. Inside the bubble, the spec- 
trum turn: over, as calculated by Lee. In the downstream region (^ 20 
min. postshock) we find over the full energy range an excellent fit of 
the differential energy spectrum dJ/dE e"V/Vp where v = 4.9 x 10® 
cm s”^ in agreement with the calculations of Forman (1981)? 

c) We measure for 35-56 keV protons an e-folding distance of 3.3 x 10^° 
cm in the shock frame (V^ij ^ 700 km s*^) in the last 10 min. upstream of 
the shock. This is in fair agreement with the value which we calculate 
from Lee's theory as L^ = 1.2 x 10*° cm using the shock (Gof^ = 40°, 

nnp = 3.1 cm ®) and particle parameters (fp(30 keV,0') = 3.3 x 10“^® cm”® 
s*, y = 1.1) of this study. The observed scale length Increases, as 
predicted, with energy. 

d) Lee (1983) predicts that upstream of the shock ion cyclotron waves 
should be generated by resonant cyclotron interactions with the energetic 
particle distribution. He estimates that the rms amplitude of the ion- 
excited waves at the shock front, normalised to that of the upstream 
field, is about 0.13. We find from ISEE-3 (6B)rms/B = 0.14 - 0.16, in 
excellent agreement with theory. 
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Abstract 


We study the quas1*para11e1 Interplanetary shock of November 11-12, 1978 
from both the collisionless shock and energetic particle points of view, using 
measurements of the Interplanetary magnetic and electric fields, solar wind 
electrons, plasma and MHD waves, and Intermediate and high-energy Ions obtained 
on ISEE-1, -2, and -3. We will characterize the Interplanetary environment 
through which the shock was propagating when It encountered the three spacecraft, 
document as completely as possible the observations of this shock, and test 
current theories of quasi-parallel shock structure and particle acceleration. 
These observations tend to confirm present self-consistent theories of first- 
order Fermi -acceleration by shocks and of col 11 si onl ess shock dissipation 
Involving firehose Instability. 10-50 keV protons accelerated by the shock 
achieved a 8 of order unity just upstream of the shock and, next to the solar 
wind energy flux, pv®/2, were the most significant dynamical parameter character- 
izing the upstream flow. The 5 keV "seed" or "upstream" protons were highly 
anisotropic and may have been accelerated by the Ion cyclotron electron heat 
flux Instability. The shock was propagating In a closed magnetic loop when 
It encountered ISEE-1, -2, and -3. 


1. INTRODUCTION 


1.1 Introductory Remarks 

The shock of November 11-12» 1978 has been characterized a*S "the 
granddaddy of all quasi-parallel Interplanetarles" (J. T. Gosling, private 
communication, 1982). Whether or not this Is apt, the shock has attracted 
considerable Interest. The facts that It was high-speed (600-700 km/sec) 
and quasi-parallel were enough to single It out for attention. It was the 
first Interplanetary shock for which an extended region of Ion acoustic 
turbulence, essentially Identical to that upstream of the bow shock, was 
found (Kennel et ^. , 1982). The region of Ion acoustic disturbance ahead 
of It proved to be surprisingly large — large enough to be of Interest to 
cosmic ray acceleration theories — and Inmedlately raised the question of 
whether other signatures of the earth's foreshock would be found equally far 
upstream of It. In the meantime, the energetic particle community was also 
looking Into the November 11-12, 1918 shock, for It was one of the strongest 
ESP particle acceleration events observed In the Intermediate energy range, 

-10 keV - 2 MeV (Reinhard , 1981). The number of papers concerned with 

the November 11-12, 1978 shock, wholly or In part. Is Increasing rapidly 
(Russell 1982; Tsurutani et^., 1983; Scholer^al^., 1983; Lee, 1983) 

and shows no signs of early abatement. 

In view of the manifest experimental and theoretical Interest In the 
November 11-12, 1978 shock, we believe a synoptic survey of the observations 
of this shock Is now In order. This paper assembles ISEE-1, -2, and -3 
observations of the Interplanetary magnetic and electric fields,' plasma, 
plasma waves, and Intermediate and high-energy Ions from 1800 UT on November 11, 


1978 to 2400 UT on November 12. 1978. a period which Includes the passage 
over the three spacecraft of a large region of upstream Interplanetary 
disturbance, the shock, multiple neutral sheet crossings downstream, and. 
finally, the piston. We will collect enough information to begin- to examine, 
for the first time, a quasi-parallel shock from the points of view of both 
Its colli si onl ess shock structure and of Its Interaction with energetic 
particles. For this shock, we will see that the two approaches cannot be 
separated. 

Our objectives In this paper are three-fold: First, to characterize 

the interplanetary environment through which the shock was propagating when 
It encountered ISEE-1. -2. and -3; second, to document as completely as 
possible the observations of the shock; and. third, using this documentation, 
to test at least the conceptual outlines and. where possible, the quantitative 
predictions of current theories of quasi-parallel shock structure and particle 
acceleration. Because of Its objectives, this paper will be long, and Its 
structure will be complex. Section 1.2 lists the Instruments whose data we 
use. Chapter 2. In which we summarize the observational and theoretical 
context from which this study has emerged. Introduces the questions upon which 
our observations have a bearing. Chapter 3 characterizes the gross properties 
of the November 11-12. 1978 shock event in terms of eight energetic particle 
and plasma diagnostics recorded during a 30-hour period containing the shock. 

In Chapter 4. we present ISEE-1 and -3 measurements of the shock and the shock 
normal solution emerging from these measurements. Chapter 4 also characterizes 
the events that occurred In the first two to three hours after shock encounter. 
Chapter 5 focuses on ISEE-3 measurements taken during the 6-hour period prior 
to shock encounter, and Chapter 6 summarizes the corresponding ISEE-1 data. 


Chapter 7 presents some experimental data bearing upon the Ion acoustic 
turbulence observed upstream of the shock and on the possible existence of 
the Ion cyclotron electron heat flux Instability upstream of the shock. 
Chapter 8 compares the theoretical predictions discussed In Section 1.2 
with our observations, and Chapter 9 closes with a few brief remarks. 


1.2 Instruments 


Table 1 lists the Instruments whose data we have assembled for this 
study. These Instruments and their capabilities have been described In a 
special Issue of the IEEE Transactions on Geoscience Electronics (Vol. GE-16, 
No. 3, July, 1978) and In subsequent scientific publications by members of 
the Individual scientific teams. The second column In Table 1 lists the 
first author and the page number In Geoscience Electronics of the article 
that describes each Instriinent. 

This paper also makes use of data on a*part1c1e count rates, obtained 
by the ISEE-3 Nuclear and Ionic Charge Distribution Particle Experiment 
(Hovestadt p. 166, 1978), and on 50 kHz electric field voltages 

obtained by the ISEE-3 3-D1mens1ona1 Radio Mapping Experiment (Knoll et 
p. 199, 1978), which were connunlcated to us privately. 

In most cases, the data assembled here are survey data presented In 
standard display formats developed by each experimental team In the course 
of their research, and the data analyses were not optimized particularly for 
the study of the November 11-12, 1978 shock. In several Instances, the 
Individual experimental teams are presently working on refined analyses of 
the data for this and other Interplanetary shock events, which they will 
publish separately. 


TABLE 1. INSTRUMENTS 


ISEE-3 


Instrument 
Solar Wind Plasma 
Low-Energy Protons 
Plasma Waves 

Vector Helium Magnetometer 


Reference 
Bame e^ a^. 
Balogh 
Scarf ^ al, 
Frandsen et 


ISEE-2 


Instrument 

Fluxgate Magnetometer 
Energetic Particles 


Reference 
Russell (p. 
Anderson et 


ISEE-1 


Instrument 

Fluxgate Magnetometer 
Energetic Particles 
Plasma Waves 

DC and Low-Frequency Electric Fields 


Reference 
Russell (p. 
Anderson et 
Gurnett et i 
Mozer et al 


(p. 160) 

.. (p. 176) 
(p. 191) 
il. (p. 195) 


239) 

(p. 213) 


239) 

al. (p. 213) 
1. (p.* 225) 
(p. 256) 


2. OBSERVATIONAL AND THEORETICAL BACKGROUND 
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2.1 Introductory Remarks 

In recent years, the space and astrophysical plasma commuMtles have 
been occupied with different aspects of three general questions: 

e Wny are quas1*para11e1 (Q||) and quasi -perpendicular (Qj^) 
colllslonless shocks so different? 

e How do colllslonless shocks accelerate particles to high 
energies? 

e How does energetic particle acceleration affect shock 
structure? 

It Is gradually becoming clear that all three questions are Interrelated, and 
we will have to consider all three to Interpret the data we have assembled for 
the Interplanetary shock event of November 11-12, 1978. We might also add a 
fourth question to the above list which data on shocks In the solar system 
alone cannot answer: To what extent can the knowledge gained from solar system 
shock studies be successfully applied to the study of cosmic ray acceleration 
by astrophysical shocks? 


2.2 Q|| and Qj^ Shock Structure 

Bow shock studies have revealed a profound difference between Qy and 
Qj^ shocks (Formlsano, 1977, Greenstadt and Fredricks, 1979). The shock normal 
angle separates the two types; Qj^ shocks have > 45**55*, and vice versa 
for Q|| shocks. Quas1*perpend1cu1ar shocks appear to be about an Ion Larmor 
radius thick (Leroy et ^. , 1982; Livesey et ^. , 1982), as measured by Jumps 
In both the magnetic field strength and plasma density. On the other hand. 

In Q|| shocks, the magnetic field undergoes a much broader and more disorderly 
transition whose spatial scale Is difficult to determine from bow«shock measure 
ments. It Is not known whether there Is a thin density Jump embedded within 
the broad region of large-scale magnetic turbulence that characterizes Q|| 
shocks. 

The facts that Q|| shocks allow significant access upstream of particles 
that have Interacted with the shock, while Qj^ shocks do not, appear to be the 
primary observational distinction between the two. As early as 1968, we knew 
that an element of sclar wlno flow could have foreknowledge of an Impending 
bow-shock crossing If It were connected magnetically to the bow-shock surface. 
The connected region of upstream disturbance has come to be known as the 
foreshock. The Interplanetary field line that Is Instantaneously tangent to 
the bow-shock surface defines the leading edge of the foreshock; at the point 
of tangr.ncy, the shock normal angle Is 90*. Near the point cf tangency, 
the locally Qj^ shock evidently accelerates electrons (Anderson, 1968, 1969; 
Anderson et al,. , 1981) and Ions (Gosling , 1978; Greenstadt ^ a],. , 1981) 
Into focused beams which escape upstream along f.^lc lines. When they are 
observed upstream, the electron and Ion beams may oe traced kinematically 
back to the shock. Although they both are generated where the *hock Is locally 


the electron beam Is encountered upstream of the Ion beam. The escaping 
electron and Ion d1str1but1<ms becoK progressively more diffuse downstream 
of their beam leading edges, on field lines that connect to a bow shock that 
Is more and more Q|| (Anderson et^., 1979; Gosling et^., 1978; 'Greens tadt 
^al., 1980). 

The superthermal particles In the foreshock generate a rich spectrum 
of magnetohydrodynamlc and plasma waves. Escaping electrons generate electron 
plasma waves (Scarf et 1971; Anderson ^ al,, 1981), low-frequency (- 1 Hz) 
whistler waves (Sentman et al 1962), and higher-frequency whistlers (Fairfield, 
1974). Ion acoustic waves are closely associated with superthermal Ions and 
electrons in the foreshock (Scarf 1971; Rodriguez and Gurnett, 1975; 

Anderson et al. , 1981). The measured Ion beam (Gary, 1%1) and diffuse Ion 
(Sentman ^ al*, 1981b) distributions are unstable to hydromagnetic waves, 
which achieve large amplitudes In the diffuse proton zone (Paschmann et 
1979; Greenstadt et al., 1980; Hoppe et al . . 1981). 

The 1n^)ress1ve clarification of foreshock phenomenology cited above has 
not Improved our fundamental understanding of the difference between Qj^ and Q|| 
shocks because of an ambiguity Inherent In the Interpretation of foreshock 
measurements. It has been argued that many, perhaps most, of the upstream 
Ions ccmie from the Ion foreshock beam (Same et al. , 1981). As such beam Ions 
propagate upstream, they (testablHze low-frequency electromagnetic waves 
which subsequently scatter and decelerate them. The decelerated Ions and 
the waves are blown downstream by the solar wind to fill the entire foreshock 
with waves and superthermal Ions. The waves are ultimately blown back Into 
the quasi-parallel zone of the shock surface, possibly accounting for the 
disordered magnetic structure of Qjj shocks. In this Interpretation, the Q|| 


C ' 


bODHShock structure we (^serve Is an artifact of the small radius of curvature 
of the bOM shock. Indeed, the foreshock has a spatial scale that equals or 
exceeds the boM shock's curvature radius. On the other hand, one can argue 
that superthermal Ions ought to escape naturally from plane Qy shocks (Edmlston 
a1. » 1982). In this case, some of the forestock phenomena we observe might 
be Inherent to Qy shock structure. Intteed, the superthermal Ions have an an 
energy density co^^)arable with that of the Interplanetary field, and, more 
significantly, the solar wind Is decelerated and deflected when It enters the 
foreshock (Bonlfazl ^ al., 1980) by an amount compatible with the momentum 
flux carried by shock-escaping Ions (Sentman et al., 1981a), suggesting that 
part of the shock transition Is accomplished in the foreshock. 

Although It Is not possible to ascertain the true extent of Qy shocks 
from bow-shock data, one Is more likely to succeed with interplanetary shocks, 
whose radii of curvature are 250-25(K) times that of the bw shock. To 
ascertain whether Qy Interplanetary shocks have foreshocks, one should begin 
by comparing measurements made at equal distances upstream of Interplanetary 
shocks and the bow shock. This in^lies searching for foreshock signatures a 
few minutes before an Interplanetary shock encounter, when the high-speed shock 
is a few earth radii from the spacecraft. At present, the search for inter- 
planetary foreshock phenomena Is Incomplete. However, Russell and Hoppe (1982) 
and Tsurutani et (1983) have found hydromagnetic waves, whose amplitudes 
and frequencies are similar to those In the earth's foreshock, a few earth 
radii ahead of Qy Interplanetary shocks. Moreover, Kennel et (1982) found 
that the an^lltude and spectrun of Ion acoustic waves a few ahead of Inter- 
planetary shocks are renarkably similar to those In the earth's foreshock. 

Thus, at least two features characteristic of the earth's foreshock also occur 


upstream of Q|| Interplanetary shocks. 

The ion acoustic waves were observed to extend several hundred earth 
radii aiwad of Q|| interplanetary shocks* the first indication that their 
foreshocks might be very much larger than is possible to infer from bow-shock 
studies (Kennel et 1982). It is of obvious Interest to Inquire whether 
other phenomena characteristic of the earth's foreshock also occur far upstream 
of interplartttary Qjj shocks. ikMvever, it would be difficult to relate a burst 
of superthermal iwis, or an outbreak of magnetic turbulence, that occurs 
several hours before an Interplanetary shock encounter* to the shock, unless 
one could also associate it with the ensend>le of phenomena characteristic of 
the earth's foreshock. Thus* case studies of interplanetary shocks are called 
for. This paper reports the results of one such study. About three hours 
before the November 11-12, 1978 quasi-parallel shock passed over ISEE-3* ion 
acoustic turbulence began to increase in intensity, at first irregularly and 
then steadily until the shock was encountered (Kennel et ^.* 1982). One 
objective of this paper will be to docmnent, with as many plasma diagnostics 
as possible* the extended region of disturbance upstream of this shock. 

Let us now briefly survey the develojxnent of our theoretical ideas 
concerning the structure of quasi-parallel shocks. It has been popular to 
separate the structure into a local shock layer* whose thickness scales as 
the Larmor radius of a thermal ion* and a more extended region, of as yet 
undetermined thickness* in which part of the shock dissipation required by 
the Rankine-Hugoniot conditions is accomplished. It is in the extended foreshock 
that the processes thought responsible for energetic particle acceleration 
occur. In other words* the overall Q|| shock consists of a foresliock and a 
plasma subshock. In general, laboratory and space plasma theoreticians have 


concentrated on the subshock and the near foreshock, and cosmic ray and 
astrophysical plasma theorists have focused on the foreshock and neglected 
subshock structure. 

Parker (1961) was the first to recognize implicitly the role of 
escaping ions in the dynamics of Q|| shocks, when he argued that parallel 
ion beams would be firehose unstable and thereby produce large-an^litude 
Alfv4n turbulence that accon^lishes the shock transition on a scale of many 
ion Larmor radii. Moiseev and Sagdeev (1963) developed a model in which 
upstream ions could be reflected from an electrostatic potential structure 
at low sonic Mach numbers; at sonic Mach numbers exceeding 1.6, the steady 
potential structure would no longer exist. Moiseev and Sagdeev (1^3) then 
argued that turbulent heating would produce a firehose unstable turbulent 
anisotropy if the upstream plasma p were sufficiently high. This point of 
view was subsequently taken up by Kennel and Sagdeev (1967) and Kennel and 
Petschek (1968), who developed a theory of low Mach number firehose shocks 
in high 0 plasmas. Auer and Volk's (1973) numerical simulation subsequently 
confirmed the general outlines of firehose shock theory. These models of the 
shock failed to recognize the in^ortance of Parker's (1961) suggestion and 
did not include the effects of escaping upstream ions. Kennel (1981) suggested 
that a fusion of the ion heat-flux and anisotropy firehose models might be 
promising. 

The above models considered only the long wavelength limit of the 
firehose instability, where it is non-resonant. On the other hand, the same 
mode is resonant for wavelengths near the thermal ion Larmor radius, or when 
the plasma 0 is less than unity (Kennel and Scarf, 1968). Recent work on 
resonant anisotropy instability has focused on the foreshock and not on the 


subshock* Gary (1961) and Sentman rt a1 . (1^1) showed that the ions escaping 
from the terrestrial subshock are unstable to the resonant Instability and, If 
the conditions are appropriate, to the non-resonant instability as well. It 
Is generally believed that this instability Is responsible for the large-amplitude 
upstream waves In the earth's foreshock. Lee (1982) has developed a self- 
consistent theory for the decay of an Ion beam escaping frcm a parallel subshock 
due to resonant quasl-l Inear scattering by low-frequency electromagnetic waves 
and the subsequent energization of Ions by scattering and shock c(»npress1on. 

In this theory, the physics of the escaping superthermal Ions and cosmic ray 
acceleration are both consequences of turbulence generated by pitch-angle 
anisotropies. 

Until recently, there vmis no evidence for waves well upstream of 
interplanetary shocks. However, Russell et al,. (1982) and Tsurutani ^ a^. (1983) 
have recently found them upstream of quasi-parallel shocks. Lee (1983) extended 
his self-consistent theory to the Interplanetary case and applied it to the 
Noveirt)er 11-12, 1978 shock. Starting with Scholar ^^'s (1983) measured 
30 keV/Q ion intensity, he was able to account for their particle measurements 
at higher energy ard to predict a wave an^jlitude and spectrum that we will 
C(»npare with observation in this paper. 

Numerical simulations have contributed substantially to quasi -parallel 
shock theory. Because of the limitations of spatial scale, numerical simulations 
treat only the sutshock. Biskamp and Weller (1972) proposed an electrostatic, 
rather than electromagnetic, ion beam instability as the dissipation mechanism 
for the strong quasi-parallel shock they simulated. Forslund and Friedberg 
(1971) and Forslund £t al^. (1980) found that the Q|j shock produces a large- 
amplitude standing whistler which decays into an ion acoustic wave and another 


whistler, and proposed that the ion wave accomplishes the Ion heating In 
the quasi -parallel shock, itecent 2-D simulations (Quest et al.» 1983) confirm 
that large-an^lltude whistler turbulence on the Ion Inertial scale length Is 
^nerated In the shock. They also find that Intense fluxes of Ions are 
reflected upstream. Kan and Swift (1983) have simulated Qy shocks In one 
spatial dimension. They find that a tdilstler wave train standing upstream 
of the shock resonantly scatters the Inccxnlng thermal Ions, and that non- 
resonant firehose tirades are created downstream. 

In sunmary, nearly all theories of Qy shock structure agree that 
large-an^lltude magnetic turbulence, with frequencies that span the range 
from well below to somewhat above the Ion cyclotron frequency. Is central 
to the dissipation In the pla^a subshock and to the dynamics of the foreshock 
ahead of It. 


2.3 Shock Acceleration of bierqetic Particles 


(Mtll recently* nrjst theories of cosmic rey acceleration have 
concentrated on elucidating how single particles can attain high energy 
by singie or multiple encounters with colllslonless shocks which are 
considered to be Infinitely thin amt whose plasma structure Is therefore 
assinned to be relatively unimportant. Looked at In this fashion* shocks 
can accelerate particles In several ways. Ions whose Larmor radius exceeds 
the shock thickness conserve their gyrophase averaged magnetic iiK)ment 
(E. N. Parker* unpublished manuscript* 1958; Chen and Armstrong* 1972; 
Shebanskll* 1972; Passes* 1979; Teresawa* 1979a*b). Such Ions approaching 
the shock from upstream would therefore be either reflected from or transmitted 
through the jump in magnetic field at the shock* depending upon their pitch 
angle. Reflected ions grad-B drfit parallel to the flow electric field* and 
they thereby acquire energy nK>re efficiently the more quasi-perpendicular 
the shock. However* since multiple reflections are needed to account for 
the observed acceleration by interplanetary shocks (Pesses* 1979)* reflected 
ions would have to scatter from upstream HHD turbulence back towards the 
shock. They then can be either re- reflected or retransmitted at their next 
encounter with the shock. Re-reflected particles can repeat the above cycle* 
and some can reach high energy. 

Energetic particles that are transmitted through the shock can be 
scattered by downstream magnetic turbulence back towards the shock. Such 
particles are subject to first-order Fermi -acceleration by multiple reflections 
between upstream and downstream waves that convect approximately with the 
local flow speed. The shock then serves primarily to decelerate the flow 
so that the scattering centers appear to converge toward one another 


In the shock frame. In the test particle limit, this mechanism does not take 
Intc account the momentum transfer between cosmic rays and the plasma. The 
Integral spectrum calculated for non*re1at1v1st1c particles Fermi -accelerated 
by Infinite plane shocks depends only upon the ratio of upstream 'and downstream 
flow speeds (Krimsky, 1977; Axford et a1 . , 1977; Bell, 1978a,b; Blandford and 
Ostrlker, 1978). Because the particle spectral Index calculated for strong 
shocks Is dose to the observed galactic cosmic ray Index, supernova shocks 
are promising candidates to accelerate galactic cosmic rays. 

For the solar system, the theory of first-order Fermi -acceleration 
has been applied to the diffuse Ions upstream of the bow shock (Elchler, 1981; 
Teresawa, 1981; Lee^^., 1981; Forman, 1981; Ellison, 1981; Lee, 1982), 
and to the so-called ESP events associated with Interplanetary shocks, 1n which 
energetic Ions are observed to increase well before the shock encounter 
(Scholer and Morfill, 1975; Scholar 1983; Lee, 1983). Lee's (1982) 

theory predicts the energy spectra of different species reported by Ipavich 
^ d_., (1981a, b) and the rough spectrum and amplitude of the low-frequency 
waves observed upstream of the bow shock by Hoppe ^ aj[. (1981) and others. 

The observed spectrum of bow-shock diffuse particles cuts off above about 
100 keV, which may be explained by the finite extent of the bow shock. Either 
a given magnetic field line remains connected to the region where the bow shock 
Is strong for a finite time, or the particles diffuse across the magnetic 
field onto field lines which no longer Interact with the shock. Either effect 
limits the nwnber of shock crossings a particle can have and, therefore, the 
energy to which it can be accelerated. 

The field line connection time Is much larger for Interplanetary shocks 
than for the bow shock, so the first-order Fermi mechanism will have longer 


to operate. According to the simplest theory, the energetic Ion fluxes 
should Increase exponentially approaching the shock, maximize at the shock, 
and hold approximately constant downstream — features characteristic of 
ESP events. The accelerated Ions upstream should be essentially Tsotropic 
In the shook frame and Isotropic In the solar wind frame. Until very recently, 
there have been relatively few measurements of moderate energy Ions In ESP 
events In the energy range (lO's of keV) that bridges the low-energy plasma 
and “seed" particles (see below) and high-energy cosmic rays. A recent study 
of 30-150 keV/Q protons and alphas In three ESP events (Scholer et ^. , 1983} 
finds that the observed particle energy and angular distributions and spatial 
profiles are consistent with first-order Fermi -acceleration theory. One of 
Scholer ^ £l^'s (1983) three events Is the shock of November 11-12, 1978, 
and we will discuss Schol‘.'r et al 's results for this shock In due course. 


2.4 Relation of Shock Structure and Particle Acceleration 


Most theories of particle acceleration by shocks treat the cosmic rays 
as test particles. E. N. Parker was one of the first to realize that cosmic 
rays had sufficient energy density to contribute to shock structure. Uentzel 
(1971) and Axford^ai. (1977) Included cosmic rays as a second fluid In 
the shock Rankine-Hugonlot relatiwis. Cosmic rays and plasma were assumed 
to have different dissipation mechanisms* whose scale length is much larger 
for the cosmic rays than for the plasma. In these models* energetic particles 
are presumed to reside upstream prior to shock passage; the cosmic ray flux* 
a free parameter* Is conserved across the shock* and if It Is large enough* 

It can create a shock entirely In cosmic rays without a subshock In the thermal 
plasma. Such multi pie- fluid shock models have been extended by Drury and Vo1k 
(1981) and Axford et al* (1982)* who modeled the plasma hydrodynamically. It 
appears that Interplanetary shocks should have both a cosmic ray foreshock 
and a plasma subshock for the range of sonic Mach numbers and particle Intensities 
encountered In the solar system. An extension of these theories to magnetohydro- 
dynamics (Hada and Kennel, 1983) confirms this conclusion. The shocks Identified 
In the Interplanetary medium should thus be considered as subshocks* and part 
of the shock transition should be accomplished across a very much more extended 
cosmic ray foreshock* an effect already observed In the terrestrial foreshock. 
Since the November 11-12* 1978 shock did have a pre-existing upstream energetic 
particle flux that was contained In a closed magnetic field region* It will be 
used to test the quantitative predictions of the above two-fluid models of shocks. 


2.5 Toward a Unified View of Shock Structure and F1rst«0rder Fermi- 
Acceleration 

Three more Interrelated questions must be answered before we arrive at 
a theory that, among other things, computes the cosmic ray Intens'lty and 
spectrum as a function of shock parameters: 

e Under what circianstances Is the Fermi *acce1 era t1 on scenario 
realized In practice? 

e How do particles that ere originally part of the thermal 
plasma reach the energy threshold where Fermi -acceleration 
begins to operate? 

e To what extent Is the turbulence that scatters thermal plasma 
In the subshock and energetic plasma In the foreshock similar 
In character? 

Recent research has made the answer to the first of the above questions 
clear. Only quasi-parallel shocks have extended regions of MHD turbulence 
both upstream and downstream, and this fact Is related to the ease with which 
not only energetic cosmic rays, but also particles on the tall of the thermal 
distribution, can free stream across the shock. 

An answer to the second question, the so-called "seed particle problem'*, 
Is urgently needed because It will ultimately tell us the efficiency with 
which shocks can convert flow energy Into cosmic ray energy. One should not 
rely on a pre-existing flux of cosmic rays upstream, and It Is generally 
believed that at least quasi -parallel shocks self-consistently generate their 
own characteristic spectrum of cosmic rays. Elchler (1979) and Ellison (1981) 
have developed a model In which both thermal and energetic Ions Interact with 
waves In essentially the same way to produce a spatial diffusion length that 
Is proportional to the Ion Lamior radius. Such a diffusion model Is capable 


of producing • high-tntrgy tall with striking tfflcltncy. Observatlonally* 

It 1$ clear that the seed particles mrt once thermal Ions that have Interacted 
once with the shock on their way to participating In the Fermi process. In 
the case of the bow shock, these are the few keV, so-called "upstream** or 
"superthermaV* Ions that were reflected from (Bame et al., 1980) or transmitted 
through (Edmiston et 1982) the shock. There have been few studies of 
the few keV analog of upstream particles from Interplanetary shocks. 

Nearly all theories of seed particles, cosmic ray acceleration, and, 
to a lesser extent, shock structure overlook the possible role of electrons 
and electron heat fluxes, not only In regulating the potential drop across 
the shock, but In participating In Instabilities that subsequently energize 
or scatter Ions. We will present evidence that an Important Interaction with 
the electrwi heat flux occurred In the foreshock far upstream of the shock of 
November 11-12, 1978. 

Existing research suggests that It may be possible to accomplish the 
entire qua si -parallel foreshock-subshock transition with low-frequency electro- 
magnetic waves as the dominant scattering mechanism. Ellison's (1^1) Monte 
Carlo simulation was based upon this presunptlon. The theories cited above 
for the resonant scattering of upstream particles (Lee, \96Z, 1983) and the 
tl^rmallzatlon of the plasma (Parker, 1961; Kenr.el and Sagdeev, 1967; Kennel 
and Petschek, 1968; Auer and Volk, 1973) all Invoke long-wavelength electro- 
magnetic waves destabilized by pitch-angle anisotropy. Even the thermallzatlw 
of Incoming plasma Ions by standing whistlers fotaid In numerical simulations 
(Que^t et al., 1983; Kan and Swift, 1963) Is a variation on the same theme. 
However, the role. If any, of short-wavelength electrostatic turbulence In 
quasi-parallel shock structure has not been defined experimentally or considered 
theoretically with any cogency. 


3. INTRODUCTORY VIEW OF NOVEMBER H-12. 1978 SHOCK 


Figurt 1 presents selected ISEE-3 date for the period 1800 UT on 
November 11, 1978 to 2400 UT on Novesber 12» 1978. The shock encountered 
ISEE-3 et 0028 UT on November 12, 1978. The data shown are, starting at the 
top panel and working downwards, 30-36 keV/Q a-part1c1e cou^t rates (M. Scholer, 
private communication, 1982), 35-56 keV proton Intensities, 50 kHz electric 
field voltages (J. L. Steinberg, private comnutilcatlon, 1982), 3 kHz electric 
field spectral ans)11tudes, 17.8 Hz magnetic field spectral amplitudes, the 
solar wind speed and density (Bane et al,., 1981), and the Interplanetary 
magnetic field magnitude. 

During the entire period from 1800 UT until the sN>ck, the Interplanetary 
field magnitude and solar wind speed remained virtually constant at 7y and 
370 km-s*', respectively. Variable fluxes of energetic protons streaming 
along magnetic field lines away from the earth Indicate that ISEE-3 had been 
magnetically connected to the earth's bow shock prior to 1725 UT on Noven^er 11. 
During the variable proton event and until 1840 UT, the superthermal electron 
fluxes were also highly variable, with their heat flux alternately toward and 
away from earth, suggesting that the period 1725-1840 UT vms one of ccm^lex, 
variable connection to the bow shock. We will concentrate on the more orcterly 
data after 1840 UT, where we will argue there was no evidence of a bow* shock 
connection. 

The solar wind density maximized, and the 35-56 keV proton Intensity 
and 3.16 kHz electric field amplitude minimized, at about 1920 UT. Thereafter, 
the 35-56 keV proton Intensity and the 3.16 kHz electric field an^lltudes 
Increased, and the solar wind density decreased, all three with superposed 
fine structure, until the shock. The five-hour, more or less steady Increase 


In the energetic proton Intensity prior to the shock Identifies this as an 
ESP event. 

It is particularly Important to note the sudden drop In solar wind 
density at 2230 UT, which Is co1nc1d«fit with sudden Increases to*p1ateaus In 
the 35*56 keV proton and 3.16 kHz Ion acoustic wave Intensities. According 
to Kennel (1^)* the electric field amplitudes between 1 and 10 kHz 

jumped near 2230 UT. Note that the 50 kHz electric field voltage diminished 
slightly at 2230 UT. The count rate of 30*36 keV/Q Helium Ions reached a 
small secondary maximum after 2230 UT. A number of other measurements, to 
be presented later, will Indicate that the state of the solar wind after 
2230 UT was distinctly different from the one preceding 2230 UT. 

The 35*56 keV proton Intensity Increased by a factor 10 In the last 
half*hour before the shock passed over ISEE*3. If this energetic proton front 
co*moved with the shock, which had a speed of about 700 km-s"^, its exponential 
scale length was about 50 R^. The measured 35*56 keV proton Intensity was 
about 3x10** p*cm'^*s"^*sr"^ just prior to the shock encounter. The corresponding 
35*56 keV proton number and energies were about 10"* p*cm"* and 0.7 ergs*cm"*, 
respectively. This energy density is probably an underestimate, because the 
measured spectrum was quite hard, but even It Is about 1/3 the energy density 
of the magnetic field upstream, 2x10"*® ergs*cm"*. Extrapolating the observed 
E“* differential proton spectrimi from 35 keV to 10 keV would Increase the 
energetic proton energy density by a factor 3.6. The fact that the energetic 
proton energy density was comparable to that of the plasma and magnetic field 
upstream of the shock suggests that energetic protons played a significant role 
in the overall structure of the shock. 

The 3.16 kHz ion acoustic wave amplitudes maximized at the shock and then 


dropped suddenly. The an^lltude of 17.8 Hz whistler mode magnetic turbulence* 
which was near threshold upstream* Increased suddenly and gradually returned 
to threshold 18 hours later. The solitudes of 0.1 - 1 kHz electric fields 
(not shown)* which were also near threshold upstream, jiai^d sudtenly at the 
shock and persisted foi* many hours dcnumstream (Kennel ^al,. , 1982). The 
50 kHz electric field voltage Jim|>ed at the shock and only returned to the 
presN>ck level after 1600 UT. Persistent enharKements of broadband low- 
frequency turbulence, 0.1 - 1 kHz electric field turbulence, and electric 
field noise above the plasma frequency (some tens of kHz In the downstream 
region), are characteristic features of the flow downstream of interplanetary 
shocks (Kennel et^., 1^2; Hoang et al . , J980). Grigorieva and Slysh (1970), 
Meyer-Vemet (1979), and Hoang et al . (1980) interpret the high-frequency 
electric field continuum as due to enhanced thermal electron plasira wave 
fluctuations. 

The count rate of 30-36 keV/Q Helium ions* which had hovered slightly 
above threshold at the shock, jt^ped by a factor 40 at the shock and gradually 
returned to threshold over the next 18 hours. The 35-56 keV proton intensity 
rose to 10* em’^-s’^-sr’^-keV'* at the shock. Note the secondary peak, which 
Is the absolute maximum, about two hours behind the shock. This peak appears 
to be associated with a local mlnlmtMn In the post-shock solar wind speed and, 
possibly, with tvx> magnetic neutral sheets which passed over ISEE-3 about when 
the 35-56 proton Intensity maximized. After the maximum* the 35-56 keV proton 
flux renalned more or less constant until about 1300 UT on November 12* at 
which time It commenced a more rapid decline and dropped suddenly near 1715 UT. 

The flare ejecta arrived at about 1715 UT* at which time the solar 
wind proton tenperature decreased, and the electron temperature diminished to 


4x10"® K, and the thermal plasma was enriched in o-particles (Bame et 1981). 
This general behavior persisted until about 0200 UT the next day. November 13. 
1978. Strong bi-directional electron streaming observed at energies exceeding 
80 eV was coincident with the arrival of the flare ejecta at ISEE-3 (Bame et al . , 
1981). The energetic protons. E > 35 keV. also exhibited bi-directional 
streaming at this time (Sanderson al^. , 1981). Bame^^. (1981) argued 
that their data were best explained if the local interplanetary magnetic field 
were part of either a magnetic bottle rooted at the sun or a disconnected loop 
propagating outward. We will suggest that the shock propagated over several 
magnetic bubbles before it arrived at ISEE-3. since several magnetic neutral 
regions were detected in the post-shock flow, and was in another bubble when 
it encountered ISEE-3. The low time resolution interplanetary field data shown 
in Figure 1 has a number of depressions in magnetic field strength downstream 
of the shock, many of which, higher resolution data will indicate, were neutral 
sheet crossings. 


4. SHOCK GEOMETRY AND DOWNSTREAM MAGNETIC STRUCTURE 


4.1 Upstream Interplanetary Field 

This paper will concentrate on the period 1800-0028 UT, when ISEE-3 
was upstream of the onc(»n1ng shock, the 35-56 keV proton Intensity Increased 
by a factor 200 , and, as we shall see, there was considerable activity In 
solar wind electrons. Ion acoustic waves, and low-frequency magnetohydrodynamic 
waves. The average properties of the Interplanetary magnetic field (IMF), 
hcwever, remained relatively constant. The magnetic field strength at ISEE-3 
was about 7v from 1800 UT on November 11, and was directed about 20*’ south 
of the ecliptic after 1840, until the shock. The IMF was approximately 
contained In the ( 6 ^, 82 ) plane that Is perpendicular to the ecliptic and 
contains the sun. The average IMF measured at ISEE-1 had roughly the same 
magnitude, and was approximately parallel to that at ISEE-3 from 1820 UT, 
when ISEE-1 exited the earth's bow shock, until Its 0058 UT interplanetary 
shock encounter. 

The stability of the IMF prior to shock encounter provided Greenstadt 
et (1982) with a good opportunity to reconstruct the magnetic field lines 
near ISEE-3 from a time-series of Individual vector measurements for the period 
after 2200 UT until the Interplanetary shock encounter. They concluded that 
the field lines at ISEE-3 passed about 60 R^ north of the ecliptic at earth 
during the entire time period, suggesting that ISEE-3 was not magnetically 
connected to the earth's bow shock. We note that the IMF lacked a significant 
"gardenhose" component throughout the Interval after 1840 UT that ISEE-3 was 
able tc observe the region upstream of the Interplanetary shock without 
confusion from effects associated with the earth's bow shock. 


4.2 Shock Encounter 


Figure 2 shows eight minutes of high>reso1ution magnetic field data 
surrounding the ISEE-3 (top panel) and ISEE-1 (bottom panel) shock encounters 
at 0028 and 0058 UT, respectively. Tte magnitude of the magnetic field, |B|, 
jumps at the shock by the same factor of 2.2 at each spacecraft. The lSEE-3 
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upstream and downstream plasma number densities are n^ s 3 cm , n 2 s 10 cm , 
and the corresponding flow speeds and proton temperatures are Vi = 375 km-s“^, 

V 2 = 575 km-s“^, and Tp^ = 3x10**® K, Tp^ = 10*® K, respectively (Bame et al^., 

1981). From the ISEE-3 data presented in this paper, we find the magnetic field 
strengths are Bi * 7.5y, 82 * 16. 5y; the magnetic field conponents will be presented 
in Figure 5. The average electron temperatures, summing over core and halo 
components, were Tg^ = 2.5x10*® K and T^^ = 4.5x10*® K. Since halo electrons 
may free stream through the shock, it may be more appropriate to consider the 
core electron temperatures, which were T^^ » 1.2x10*® K and T^^ a 2.5x10*® K. 

Ti« solar wind accelerated downstream of the shock and reached a peak speed 
approaching 800 km-s’^ one-half to one hour after the shock encounter. 

If this interplanetary shock originated in a 2N flare observed near the 
central meridian of t^w sun at 0100 UT on November 10, 1978, the average shock 
velocity was 880 km-s’^ (Reinhard ^ a^. , 1981). The 30-m1nute delay between 
the ISEE-3 and ISEE -1 shock encounters indicates that the shock speed was about 
700 km-s"^ near the earth. 

We have calculated shock normals using four different techniques: 

( 1 ) Using coplanarity across the shock at ISEE-1; (2) using coplanarity at 
ISEE-3; (3) using three-dimensional solar wind data in the mixed mode method 
of Abraham-Schrauner and Yun (1976); and (4) taking the cross-product of the 
field jumps observed at ISEE-1 and -3 in various con4>inations (the double 


(telta B technique). The results using method (3) using 2-d1mens1onal data, 
which were reported by Kennel et iH. (1982), will be amended here. The shock 
solutions from these four techniques are presented In Table 2. The first 
column In Table 2 shows the components of the shock normal n In GSM coordinates; 
the second column shows the components of the magnetic field vector B; and 
the third colimn shows the shock normal * cos“‘ n*B/lBl. The double delta 
B row contains thr^e entries In the last two colianns, corresponding to various 
estimates of the upstream magnetic field at ISEE-1 and -3. 

All four techniques give similar shock normals and Indicate that 
was roughly 40^. Using method (3), we have determined that the component of 
the shock velocity normal to the shock was 609 km-s”^. For the upstream 
parameters In the first paragraph of this section, the ratio of the 
electron pressure, based upon the average electron temperature, to the magnetic 
pressure was 0.54. The thermal proton was about 0.065. For a total upstream 
Be Bp of 0.55, solutiwi of the Rankine-Hugonlot relations Indicates that a 
fast Mach number Mp of 2.6 Is necessary to give a jump In field strength of 
3.2 In the shock plane and a jump of 2.2 In the total field, approximately 
what was observed. Thus, the interplanetary shock of Noven^er 12, 1978 was 
of moderate strength, with 6g|^ near the Qjj - Qj^ transition (Greenstadt and 
Fredricks, 1979). 

The shock normal and Rankine-Hugonlot solutions In this paper and, to 
our knowledge, all others In the literature to date do not take Into account 
the effects of energetic Ions, whose B was close to unity and possibly exceeded 
that of the plasma just upstream of the November 12, 1978 shock. These solutions 
are, strictly speaking, for the local subshock. The full Rankine-Hugonlot 
conditions should be taken between the points about a half-hour upstream where 
the energetic protons began to Increase rapidly In Intensity and a point well 
downstream of the subshock. 


Figures 3a and 3b sketch the geometry of the shock at the moment It 
passed over ISEE-3 at 0028 UT on November 12, 1978. We used the method (3) 
shock normal solution discussed above and lSEE-1 and -3 Interplanetary field 
data. Taken by Itself, the X-Y projection (Figure 3b) might suggest that both 
ISEE-1 and ISEE>3 could have been ccmnected simultaneously to the Interplanetary 
shock and the bow shock, because the By-coit^onent was nearly zero. However, 
the X-Z projection (Figure 3a) suggests that the IMF had a southward component 
that would cause the ISEE-3 field line to miss the bow shock In the absence 
of significant field line curvature between ISEE-3 and the bow shock, as 
Greenstadt et al . (1982) argued. Tlie fact that the simultaneously-measured 
ISEE-1 and -3 magnetic fields were B»re or less parallel suggests that such 
curvature was indeed absent at the time of ISEE-3 encounter. The fields at 
ISEE-1 and -3 were approximately parallel from at least 2100 UT on November 11, 

1978 until the shock was encountered at each spacecraft. A time projection 
of the ISEE-3 magnetic field line indicates that It missed the bow shock throughout 
this period (Greenstadt et al., 1982); however. It Is likely that ISEE-1 was 
connected to both shocks. The uniformity of the IMF indicates that the structure 
In the Interplanetary field had a spatial scale exceeding the 200 R£ that separate 
ISEE-1 and ISEE-3. 


TABLE 2 


Technique 

Nomal 


Uostream B(y) 

^BN 


Hjl Hi 

!z 

h ^ h 


ISEE-1 Co-planarity 

-.966, .W3, 

-.256 

6.13, 0.31, -3.05 

41.6 

ISEE-3 Co-planarity 

-.942, .298, 

-.152 

5.97, -1.07, -2.82 

34.9 

ISEE-3 Mixed Mode 

-.966, .073, 

-.247 

6.40. -0.99, -3.23 

41.0 

Double Delta B 

-.954, .104, 

-.279 

6.13, 0.31, -3.05 

43.5' 




5.97, -1.07, -2.82 

41.3' 




6.40, -0.99, -3.23 

42.8' 


4.3 Properties of MHD Waves Near the Shock 


Tsurutani ^ (1983) have studied the polarizations of the magnetic 

field oscillations measured by ISEE-3 In assxiatlon with the Nov^er 11-12, 
1978 shock event. We summarize their conclusions here. Large-amplitude 
compresslonal waves with about 7-second periods are apparent just ahead of 
the shock In the top panel of Figure 2. The waves between 0026:20 and 0027 UT 
were found to be elllptically polarized. A minimum variance analysis showed 
that the angle 6|^ between the wave propagation vector t and the local average 
magnetic field direction § was 7®. The By-component in the top panel of 
Figure 2 revealed shorter-period (3 sec) waves of even larger amplitude just 
downstream of the shock. The waves measured between 0028:19 and 0028:30 UT 
were right-hand elllptically polarized with 6|^ * 12®. 20-second quasi-periodic 
waves occur in the interval 0029:20 - 0031 UT; these waves were non -compressive, 
predominantly right-hand elliptically polarized with 6j^ * 5®. Thus, in all 
three cases, the waves near the shock propagated essentially parallel to the 
local magnetic field. 


4.4 Rotational Discontinuity 


As Figure 2 shows, the shock Initiated a period of Intense magnetic 
oscillations that was terminated by a strong rotational discontinuity six to 
seven minutes later observed at both spacecraft. The discontinuity rotated 
the average IMF from the X-Z plane configuration that had prevailed for the 
previous six hours to an "antl-gardenhose" configuration In the X-Y plane. 

Figure 4 shows two components of the electric field and three components 
of the IMF measured at ISEE-1 between 0050 UT and 0115 UT, a period which 
contains the shock at 0058 UT and the rotational discontinuity at 0105 UT. 

These electric field components are accurate to about ± 1 mV/m, and the third 
component may be estimated assuming E*B ^ 0. kle have verified that these 
electric fields are consistent with the solar wind speed measured upstream 
of the shock. The general anti correlation of Ey and 82 Is consistent with 
an MHO flow In the -X direction, and a strong rotation of the electric field 
Is apparent at the magnetic discontinuity. Note the partial re-entry Into the 
state upstream of the rotational discontinuity at about 0108 UT. 

The large-amplitude magnetic oscillations In Figure 2, which Figure 4 
shows were essentially linearly polarized, were concentrated between the shock 
and the rotational discontinuity. 


4.5 Downstream Neutral Region? 


ISEE-3 and ISEE>1 each encountered a succession of magnetic neutral 
sheets downstream of the shock. Figure 5 shows the first of these, which 
followed 26 and 24 minutes after the shock passed over lSEE-3 and -1, 
respectively. The top four panels show three components and the magnitude 
of the IMF iMasured at ISEE-3, and the bottom panel shows the IMF magnitude 
at ISEE-1. The IMF magnitude reached a mlnimwn of less than 1/ at about 
0054:20 at ISEE-3 and a minimum of about 2.5y at 0122:20 at ISEE-1. The 
28-m1nute delay between the neutral sheet encounters Is roughly consistent 
with the solar wind travel time between ISEE-3 and -1. The ISEE-1 neutral 
sheet encounter differed In one significant way from Its ISEE-3 counterpart. 

The second panel down In Figure 5 Indicates that at ISEE-3, the By-component 
of the IMF changed sign, whereas at ISEE-1 our data (not shown) Indicate 
that By dipped to zero and recovered back to Its original positive value 
during the neutral sheet crossing. 

Figure 6 shows the magnitude of the IMF for the next three partial 
neutral sheet encounters by ISEE-3, near 0114 UT, 0207 UT, and 0229 UT on 
November 12, 1976. 

In summary, the Interplanetary magnetic field was unusually configured 
during the entire Noven^er 11-12, 1978 shock event. From 1800 UT on November 11, 
six and one-half hours upstream of the shock until many hours downstream, the 
IMF acted as though It were not connected to the sun. The IMF was In the X-Z 
plane throughout the upstream period; a rotational discontinuity that followed 
the shock did create a non-zero By, but it rotated the IMF Into an antl-gardenhose 
sense, with both B^^ and By positive. These observations may have a simple 
explanation, for If ISEE-1 and -3 ¥«re near the general field reversal of the 
solar wind, one might expect neutral sheet dissipation to disconnect the IMF 


from the sun, and there would be no reason for the disconnected field lines 
to be In the gardenhose configuration. The succession of neutral sheets that 
passed over ISEE-1 and -3 following the shock support this view and have an 
Interesting Implication. ISEE«3 and -1 may have detected a train of X and 
0>type neutral regions In the field reversal layer. The magnetic field 
surrounding an 0*type neutral region will be essentially closed. Bame et al. 
(1981) associated one such closed region with the driver of the November 12 » 1978 
shock. The shock had evidently passed over several others before aniving at 
ISEE-3, as the succession of neutral sheets suggests. Furthermore, we will 
argue here that the shock was In the process of devouring another closed magnetic 
loop, which ISEE-3 and *1 entered two hours and one hour before their shock 
encounters, respectively. 


S. Ut»STR£AH REGION: lSEE-3 


Figure 7 shows model electron perimeters derived from meesurements made 
by the LANL/MPI solar wind ulactron detector for the period 1800 UT on November 11 
to 0100 on November 1?. These model parameters are defined, and their significance 
discussed, by Feldman et (1982). The top panel of Figure 7 shows the electron 
heat flux parallel to the magnetic field, the bulk of which Is carried by ''halo'' 
electrons with energies of tens of eV. With the exception of a short Interval 
near 1830 UT, this heat flux was directed away from the sun. The middle panel 
shows electron (tensities for the same time period; the upper curve shows the 
total density, and the bottom curve shows the density of the superthermal halo 
population. The bottom panel shows mo(tel electron temperatures, with the halo 
temperatures at the top, the core at the bottom, and the density-weighted average 
ten^erature Inbetween. The temperatures have been shaded to Indicate the thermal 
anisotropy; the top of each shaded graph Is the parallel temperature, and the 
bottcxn Is the perpendicular. 

ISEE-3 was connected to the earth's bow shock for several hours prior 
to 1840 UT. The last portion of this period of connect1(m Is evident In the 
top panel of Figure 7, which shows that the heat flux fluctuated violently. 

The electron diagnostics became much more regular after 1840 UT and never fluctuated 
so violently again until the 0028 UT shock encounter. The total density (top 
curve, middle panel) rose to a local maximum of about 15 cm*’ at 1930 UT, where 
Figure 1 Indicates the 35-56 keV proton flux minimized, and settled to a plateau. 

The density dropped suddenly at 2100 UT, recovered somewhat, and then resumed 
Its decline, which ended In a sudden drop at about 2220 UT. The..dens1ty was 
constant from then until the shock encoimter. The halo electron parallel 
temperature rose steadily from 2130 to 2230 UT and held more or less constant 


until the shock. The core electron temperature rose gradually and then more 
rapidly between 2130 and 2230 UT; It, too, held constant from about 2230 UT 
until the shock encMinter. The electrcm heat flux, which had been more or 
less constant since 1900 UT, diminished between 2200 and :230 UT,* and then 
remained constant until the shxk. 

The solar wind electron parameters displayed In Figure 7 all Indicate 
that the solar wind entered a new state after 2230, In which the electron density 
and net heat flux were lower, and the electron temperature was higher, than 
before 2200 UT. Figure 8 displays refined halo electron parameters that help 
to diagnose the transition between these two states. The top panel shows the 
halo drift speed parallel to the magnetic field — the first parallel velocity 
moment of the entire halo distribution. The halo drift speed was between 
1500*2000 km/ sec, directed away from the sun, between 1840 and 2100 UTi It 
dropped to less than 1000 km/ sec and held rouchly constant between 2230 UT 
and the shock at 0028 UT. To help define the source of the halo electrons, 
we consider the properties of those electrons strMmIng along the magnetic 
field toward and away from the sun separately. The partial density and 
parallel speed moments of the "away" and "toward" semi-distribution functions 
are displayed In the middle and bottom panels of Figure 8. The density 
of "away" electrons (top curve of middle paiwl) did not decrease, but 
Increased somewhat, as the new solar wind region moved over ISEE-3. In 
addition, the "away" drift speed (t(H) curve of bottom panel) also Increased, 
so that the thermal flux carried by electrons streaming away from the sun was 
about a factor 2 larger after 2230 UT than before. On the other hand, the partial 
density and speed of the electrons streaming toward the sun (bottom curves of 
the middle and bottom panels), which were small prior to 2200 UT, both rose 
abruptly beginning about 2200 UT to values more nearly comparable with those of 


the away electrons which renained constant between 2230 UT and the shock. The 
mt heat flux and wean halo drift (top panels of Figures 7 and 8, respectively) 
did not decrease tecause the halo electrons streaming avmy from the sun and the 
micomlng Interplanetary shock diminished In Intensity, but because a new con^onent 
of electrons streaming toward the sun was present after 2230 UT, and the halo 
electron distrloutlon was more symmetrical. 

Figure 9 diagnoses the transition between the two solar wind states 
1(tent1f1ed above with higher time resolution. Its top panel shows the 35-56 keV 
protOT flux between 2200 and 23CK) UT; the middle panel displays the first parallel 
velocity moment of the 35-56 keV proton angular distribution, and the bottom 
panel presents the electric field spectral density of 3.16 kHz Ion acoustic waves. 
Prior to 2220 UT, the 35-56 keV protons streaming velocity was directed away from 
the sun, and after 2230 UT, it was essentially zero, by which time the ion acoustic 
waves had increased in intensity. The 35-56 keV ions remained isotropic until 
about 2400 UT, at which time a Significant streaming anisotropy away from the sun 
developed and increased until the shock encounter. The Ion acoustic waves also 
Increased in Intensity between 2400 UT and the shock (see Figure 10). In contrast 
to the halo electrons, the 35-56 keV protms did not exhibit a bi-directional 
streaming, but tended toward isotropy between 2230 and 2400 UT. 

The fluxes of halo electrons streaming toward the sun detected after 
2230 UT can be interpreted in at least two ways. ISEE-3 could have been magnetically 
connected to the earth's bow shock after 2230 UT and thus have detected bow-shock 
electrons, or ISEE-3 could have entered a closed magnetic loop at 2230 UT that 
Intersected the Interplanetary shock at b“'th ends. In the second case, one might 
expect ISEE-3 to detect stronger fluxes streaming away from the point of nearer 
connection and weaker fluxes appearing to stream toward the sun that were generated 
at the point of further connection and followed a long path around the closed 
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field to reach ISEE-3. Although electron data alone might not distinguish 
between these two alternatives, we believe the balance of evidence favors the 
second. The data from the period before 1840 UT, when ISEE-3 was definitely 
sporadically connected to earth, differed in internal character from that after 
2230 UT. The electron heat flux parameters fluctuated violently, and the heat 
flux intermittently reversed sign, prior to 1840 UT; the > 35 keV protons also 
had a streaming anisotropy away from earth prior to 1725 UT. On the other hand, 
the electron parameters were stable, and a sunward proton anisotropy was absent, 
after 2230 UT. We will present further argynents supporting the magnetic loop 
hypothesis later. 

Figure 10 shows four hours of 3.16 kHz electric field data (top panel) 
and the rms average (.01-3 Hz) interplanetary field fluctuation amplitude 
surrounding the ISEE-3 shock encounter (bottom two panels). Shown are the 
Y and Z-components of the rms fluctuation amplitude; we have shaded the periods 
when the squared rms deviation exceeded 0.1 to guide the eye. The 3.16 kHz 
electric field and the low-frequency magnetic field fluctuations appear to track 
one another; bursts of ion acoustic noise near 2140 UT and 2200 UT were accompanied 
by bursts of hydromagnetic waves. There is a definite change of character after 
about 2230 UT, which we have denoted by a vertical line connecting all three 
panels; the times of the shock, field rotation, and first neutral sheet encounter 
are indicated in similar fashion. After 2230 UT, the ion acoustic and magneto- 
hydrodynamic wave anplitudes remained intense and nearly constant until 2400 UT. 

The ion acoustic an^litudes increased between 2400 UT and the shock encounter, 
during whi" 'i time the 35-56 keV protons increased in intensity and again had a 
significant streaming anisotropy away from the oncoming shock, the rms inter- 
planetary field amplitudes, at best, increased only slightly during this last 
half hour before shock encounter. 


We have nia<te a rough determination of the frequency dependencies of 
the polarization and an^>11tude of the H(€ turbulence that was observed during 
the last 16 minutes prior to shock encounter. Figure 11 presents our results 
for the period (X}12-0028 UT (solid line), and also for the shorter period 
0025:30-0027 UT, Immediately upstream of the shock. A11 the data have 
binned Into factor 2 Intervals of spacecraft frame frequency. Starting from 
the top, the panels In Figure 11 present the frequency dependencies of (1) the 
angle 6|^ between the wave propagation vector, K, and the IMF, (2) the 
ell Ipti city and sense of rotation of the polarlzatlw In the spacecraft frame; 

(3) the percent polarization; and (4) the total rms amplitude and 

the rms compresslona! an^lltude <(6B^j)^>, both normalized to the Interplanetary 
field magnitude. The proton gyrofrequency based upon the measured magnetic 
field strength of - 7.5 y was roughly 0.12 Hz; the maxiimim Doppler-shift frequency 
KV/2v for KRj « 1, where Rj Is the thermal proton Lannor radius based upon the 
measured upstream proton temperature, 7pi = 3x10“® K, and flow speed, 370 km-s"^, 
was roughly 1.5 Hz. 

Figure 11 Indicates that the total rms amplitude was concentrated between 
0.01 and 0.16 Hz and reached a peak normalized amplitude of 0.16 between 0.02 
and 0.08 Hz. Thus, the MHO turbulence had frequencies below the local proton 
cyclotron frequency and, since the power diminished above 0.16 Hz, wavelengths 
that were considerably larger than the thermal Ion Lanror radius. The v«ves 
propagated at about 20® to the magnetic field between 0.01 and 0.16 Hz, and. 

In keeping with this fact, their rms compresslonal power was much smaller than 
their total rms power In this frequency band. The waves were weakly, predominantly 
left hand elliptically polarized In the spacecraft frame. 


Figure 12 compares the properties of 35-1600 keV protons (top two panels) 
with the fluctuating 62-coniponent of the .01-3 Hz Interplanetary field from 
2100 UT on Noveni^r 11 to 0100 UT on November 12* 1978. The arrivals of the 
bubble, the shock, and the rotational discontinuity at ISEE-3 are Indicated by 
the vertical lines that connect all three panels. The top panel shows proton 
Intensities for eight channels, labeled as follows: 1 (35-56 keV), 2 (56-91), 

3 (91-147), 4 (147-237), 5 (237-383), 6 (383-619), 7 (619-1000), and 8 (1000-1600 
keV). "he middle panel shows the effective energy power law spectral Index obtained 
from the ratios of the differential energy fluxes In channels 1 and 2 (bottom 
trace), 3 and 4 (second from bottom), 5 and 6 (second from top), and 7 and 8 (top). 
The proton Intensities Increased upstream of the shock at all energies. Even 
the 1-1.6 MeV proton Intensities Increased by a factor 3 during the 41 hours 
prior to shock encounter. However, only the lowest energy channels, 1, 2, and 3, 
responded significantly to the arrival of the bubble. The proton Intensities in 
channels 1, 2, and 3, which dipped about 2200 UT, rose prior to 2230 UT and 
reached a plateau Inside the bubble. They diminished somewhat near 2320 UT and 
then began their rapid rise as the shock approached. The intensification was 
most pronounced at the lowest energy. The approach of a shock which is accelerating 
protons Is evident In the spectral Index plots In the middle panel. The 
spectral Index determined from the lowest-energy channels actually turned 
negative. Indicating that the peak of the energy distribution was above 56 keV, 
at 2320, before rising to Its post-shock value of about 1.3. Scholar et al . 

(1983) have shown that the energetic proton e-folding distance was roughly 
proportional to the square root of the particle energy upstream of the 
November 12, 1978 shock, so that the negative spectral Index In Figure 12 
may signify the prior arrival of a diffusion front at energies above 56 keV. 


6. UPSTREAM REGION: ISEE-1 


Figure 13 displays the Frequency-Receiver (SFR) plasma wave 
electric field measurements made on ISEE-1 between 1800 UT on Noven^er 11, 1978 
and 0200 UT on November 12, 1978. ISEE-1 entered the solar wind at approximately 
1840 UT, at which point the intense low-frequency waves (f <2 kHz) characteristic 
of the magnetosheath diminished in Intensity. As It crossed the bow shock, ISEE-1 
began to detect 40-50 kHz electron plasma oscillations which persisted without 
significant interruption until the interplanetary shock encounter at 0058 UT. 

The frequency of the electron plasma waves decreased throughout this period, 
in keeping with the general decline of electron density found after 1930 UT 
by ISEE-3. However, the more rapid decline in density that began abruptly at 
2100 UT at ISEE-3 appears to have been delayed by an hour or more at ISEE-1. 

These plasma waves, which were probably generated by energetic electrons escaping 
from the bow shock (Anderson, 1968; Anderson , 1979; Feldman al[. , 1982), 

were not detected at ISEE-3. Thus, while ISEE-1 may well have been connected 
to the earth's bow shock prior to the interplanetary shock encounter, these 
results again suggest that ISEE-3 was not. 

The 500 Hz - 5 kHz ion acoustic amplitudes at ISEE-1 increased at about 
2300 UT, approximately two hours before shock encounter. If they had been 
generated by particles escaping from the earth's bow shock* one might 
have expected them to have been present between the bow shock exit and the 
interplanetary shock encounter, for the average IMF did not change its direction 
significantly during most of this period. In fact, ion acoustic waves were 
imich less intense between 2000 UT and 2300 UT, suggesting that the increase in 
intensity after 2300 UT was related to the approach of the interplanetary shock. 

Let us suppose that ISEE-3 entered a closed magnetic bubble at 2300 UT. 


The leading edge of the bubble •• the separatrix between "open" and "closed" 
field lines — should propagate at the solar wind speed ahead of the shock, 

* 370 km/sec. One would thus expect ISEE-1 to enter the bubble about an hour 
later than ISEE-3. Figure 14 compares the rms fluctuation amplitudes of the 
total (Bj) and Y-compor«nt (By) of the interplanetary field measured at ISEE-1 
and -3 aligned in time according to the shock encounter. The estimated times 
of arrival of the bitf)ble at ISEE-1 and ISEE-3 are Indicated by an arrw. While 
the hydromagnetic wave activity is sporadic prior to the bubble arrival at each 
spacecraft, the rms anplltudes Increased at 2230 (ISEE-3) and 2330 (ISEE-l) 
and remained enhanced until the shock encounters. 

Figure 15 presents lSEE-1 data for the period 2200-0200 UT that diagnoses 
further the arrival of the bubble at ISEE-1. The top two panels show 1 kHz and 
1.78 kHz electric field amplitudes; the middle panel shows the rms B 2 interplanetary 
field amplitude; and the bottom panels show the fluxes of 2 and 6 keV protons 
measured along the spin axis of the ISEE-1 spacecraft. Prior to 2330 UT, sporadic 
bursts of 2 and 6 keV protons, magnetohydrodynamic waves, and Ion acoustic waves 
accompany one another. The change in solar wind state at 2330 UT is strikingly 
apparent in the proton fluxes, which increase abruptly and remain enhanced — 
albeit with violent fluctuations -- until the shock encounter. 

The onset of the magnetic turbulence appears to have been rather sudden 
at the times of bubble arrival at ISEE-1 and -3. Figure 16 c(xnpares high 
resolution measurements of the interplanetary field for 12 minutes surrounding 
2230 UT at ISEE-3 (top curves in each panel) and 2325 UT at ISEE-1 (bottom 
curves in each panel). The ISEE-3 an^litude scales are to the left of each 
panel, and the ISEE-1 are to the right. A clear onset Is apparent In each 
data set; the waves appear to have 15 to 20-second periods and to be largely 
transversely polarized, similar to the hydromagnetic waves upstream of the 


L y 


earth's bow shock. Since such magnetohydrodynamlc waves propagate within 
about 100 km/ sec of the solar wind speed, the sudden onsets of wave turbulence 
suggest that the waves, and the particles responsible for their generation, 
were already present In the bubble before It moved over each spacecraft. 

Figure 15 seems to support this conclusion. 


7. ORIGIN OF UPSTREAM ELECTROSTATIC WAVES 


Figure 17 presents ISEE-1 (bottom two panels) and ISEE-2 (top two 
panels) 6 keV proton fluxes and the Z-component of the IMF for the period 
2324-23^ UT surrounding the passage of the closed ma^etlc field bubble 
over the two spacecraft. The 6 keV proton fluxes exhibited factor 100 
modulations whose time scales were comparable with the periods of the 
hydromagnetic waves. A similar modulation was observed at 2 keV. The five 
vertical lines, which have been drawn In to guide the eye, indicate that 
the € keV proton fluxes peaked when 6 ^ approached zero. The proton detectors 
are sensitive to a narrow cone of particle arrival angles centered along 
the spin axes of the two spacecraft, which In turn are very nearly perpendicular 
to the (Bjj,By) plane. Therefore, when B 2 neared zero, 90® pitch-angle protons 
were being detected. We therefore conclude that the 2 and 6 keV protons 
were highly anisotropic, with a strong peak near 90® pitch angle. Similar 
modulations of the 6 keV proton flux were detected sporadically but continually 
frcxn the entry of ISEE-1 and -2 Into the bubble near 2330 UT until the shock 
encounter at 0058 UT. 

ISEE-2 was about 4000 km upstream of ISEE-1 at the time the measurements 
displayed in Figure 17 were taken. The ISEE-1 and -2 data were offset by 
12 seconds — about the solar wind travel time between the two spacecraft, 
based on a speed of 370 km/sec — to shcm the correspcmdence between the peaks 
In proton flux and 82 at ISEE-1 and -2. These data suggest that the anisotropic 
6 keV protons and hydrcmiagnetlc waves were convected downstream at about the 
solar wind speed between the two spacecraft. 

The very strong electron heat flux In the region upstream of the shock 
(Figure 7), which Is carried by halo electrons, would lead to a strong electrical 


current unless the core electrons also had a net drift relative to the Ions 
that neutralizes the heat current. The top two panels In Figure 18 plot 
the measured model parameters and (Nq/N)V|^^, where Vp^ and Vp^ 

are the halo and core drift velocities parallel to the magnetic field, 
respectively. The extent to which the two panels are mirror- symmetric 
Indicates the degree of current neutrality. Ne note that the cold electrons 
drifted with respect to the Ions with a speed of about 200 km/sec between 
2230 UT and the shock. 

The bottom panel shows 328-sec average 1.78 kHz electric fields 
(bar graph), and the peak amplitudes (Individual points) detected during 
each 328-second Interval. The 1.78 kHz Ion acoustic waves were similar 
to those at 3.16 kHz (Figure 10) upstream of the shock. We have enough data 
to ask whether these waves were due to the Ion acoustic heat flux Instability 
(Forslund, 1970). For the measured core electron tenperature (Figure 7) and 
Bame et al^'s (1981) reported temperatures, the calculations of Kindel and 
Kennel (1971) Indicate that a core electron drift of 200 km/sec will be 
stable to Ion acoustic waves. We conclude, therefore, that the upstream 
Ion acoustic waves may have been thermal fluctuations that were enhanced when 
the core electron drift Increased. On the other hand, a drift of 200 km/sec 
Is close to the Instability threshold of electrostatic Ion cyclotron waves. 

This suggests that the highly anisotropic 6 keV protons In Figure 17 may 
have been generated by a solar wind analog of the process that accelerates 
auroral Ions perpendicular to the magnetic field to form "conic" distributions 
(Okuda and Ashour-Abdalla, 1982; Ashour-Abdalla and (Ncuda, 1982). We also 
note that the 200 km-s*^ drift also exceeds the Alfven speed, which was about 
100 km-s*^, so that the electromagnetic Ion cyclotron wave might also be unstabl 
to the electron heat flux (Forslund, 1970). 


8. OBSERVATIONAL TESTS OF THEORY 


8.1 Introductory Remarks 

In this chapter, we use the data we have assembled to test some of the 
theoretical predictions discussed In Section 1.2. In Section 8.2 we compare 
the properties of the energetic particles upstream of the November 11-12, 1978 
shock with theory. We discuss the electromagnetic waves upstream of the shock 
In Section 8.3. In Section 8.4, we briefly consider the question of the 
energetic particle effects on shock structure, and In Section 8.5 we compare 
observations and theory of the state downstream of the plasma shock. 


8.2 Energetic Particles Upstream 


The theory of first-order Fermi acceleration makes the following broad 
predictions about the behavior of Ions whose energies exceed a certain threshold 
of about 10 keV: 

• The energetic proton distribution upstream should 
Increase essentially exponentially leading up to 
the shock. The exponential scale length can be 
calculated using a self-consistent model of the 
generation of HHD waves by the energetic particles. 

e The energetic protons upstream should be essentially 
Isotropic In the frame of the shock. The anisotropy 
In the solar wind frame provides the free energy 
for wave growth. 

e The energetic proton momentum distribution should 
be a pow«r lew whose spectral Index 8 Is given 
by 0’^ • 1/3(1 - V 2 /V 1 ), where Vx and V 2 are the 
solar wind velocity components perpendicular to 
the stock plane upstream and downstream, respectively. 

e The power law should prevail between a lower limit 
momentwn, determined by shock Injection processes, 
and an upper limit, set by the condition that the 
momentum diffusion time be less than the shock 
connection time. 

The four points above have been considered In detail by Lee (1983) and 
Scholar et a1 . (1983) for the shock of November 11-12, 1978. We will paraphrase 
Lee's (1983) discussion of them so that the reader may have a more complete 
perspective on the physics of this shock. 

Scholer et a^. (1983) found that protons upstream between- 15-150 keV 
Increased exponentially with a scale length that depends upon the square root 
of the particle energy. Lee (1983) derives the following expression for the 


scale length at the shock: 
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where 0 Is the momentum spectral Index, 1s the Alfven speed, Op Is 

the proton cyclotron frequency, Mp Is the proton moss, fp(E) Is the 

omnidirectional proton momentum distribution function at the energy E, 

and V»o ■ 2.4x10* cm-s*^ Is the speed of a 30 keV proton. According to 

Scholer^til. (1983), 0 ■ 4.1 and fp(30 keV) • 3x10’*'* s* cm’*. Upstream 

of the shock, we find 6 « 7.5y and N ■ 3 cm’* (Bame et ^. , 1981), so that 

s 100 km-s’*. We choose 6g|^ » 40°, In keeping with the revised Rank1ne> 

Hugonlot solutions In Section 4.2. Substituting these values Into (8.1), 

. 10 

we obtain L^(30 keV) a 1.2x10 cm, about a factor c smaller than the 

10 

measured value, 1^(30 keV) a 2.6x10 cm (Scholar ^ a]^., 1983). 

Scholar et (1983) did find that the upstream energetic proton 
distribution was approximately Isotropic In the shock frame. We agree with 
this conclusion for the 30-m1nute period before the shock crossing to which 
Lee's (1983) theory primarily applies. The matter of the upstream and downstream 
proton spectra has been discussed by Sanderson et al_. (1981) and Scholar et a][. 
(1983), and the exact form of the distribution function may be subject to scxne 
revision. Scholar et (1983), in an atten^t to close with theory, fit the 
momentum distribution •'•ith a power law of spectral Index 0 ■ 4.1. Our measured 
compression of the components of the Interplanetary field tangential to the 
shock, 3.2, agrees with the observed density compressloi of about 3 and Implies 
Va/Vi« 1/3.2. inserting this value Into the theoretical expression above leads 
to 0 > 4.36, In excellent agreement with the observational estimate. Note that 
we used the compression ratio of the plasma subshock and not the (slightly) 
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larger ratio \;hat would apply to the full plasma shock-energetic particle 
foreshock system. 

By setting the characteristic acceleration time equal to the solar 
wind travel time from the sun, Lee (1983) finds that the power laV spectrum 
should extend up to about 200 keV/nucleon, In rough agreement with observation. 


8.3 MHO Waves Upstream 

Lee's (1983) self-consistent theory mekes the foil owing predictions 
concerning the MHO waves upstream of a quasi -parallel shock: 

e Electromagnetic Ion cyclotron waves should be 
generated by resonant cyclotron Interactions with 
the energetic particle distribution. 

e The wave spectnm should peak at a frequency that 
resonates with the lowest energy particle that 
participates In the Fermi process and fall to 
the ambient solar wind background above that 
frequency. 

e The waves upstream of Interplanetary shocks 
should be more nearly unpolarized than those 
In the earth's foreshock. 

e Given the energetic particle distribution at 
one energy at the shock, the spatial and S|MCtral 
distrlbuticn of the wave amplitude may be determined 
self-consistently. 

Lee's (1983) theory was restricted to electromagnetic waves that propagate 
parallel to the ambient magnetic field. However, Tsurutani et aj_. (1983) and 
we find that the waves near the November 11-12, 1978 shock did propagate almost 
parallel to the magnetic field. Moreover, the waves upstream of Interplanetary 
shocks show a smaller polarization than those In the earth's foreshock. 

Let us estimate the freqi^ncy of the waves that resonate with an 
energetic proton whose velocity parallel to the magnetic field Is Vj^ In the 
' solar wind frame. The cyclotron resonance condition Is approximately 
where ICy Is the parallel wave nueber. Since the phase 
velocity Is much less than the solar wind speed, the frequency f observed 



In the spacecraft frme will be bounded by 


where V Is the solar wind speed and f . Is the proton cyclotron frequency 

I 

In Hz. For B » 7.5 y as observed, f^^ * 0.12 Hz. The Iwnest-energy 
particle that can escape upstream of the stock will have a speed of 
about > 7(M) km*s*^ In the spacecraft frame. It seems reasonable 
to assume that the shock Injects particles at about twice this speed, 

1400 km-s"^, corresponding to 10 keV energy In the spacecraft frame, 
since the observed particle spatial profiles are smooth above 10 keV 
(see Scholer al- , 1^3, and our Figure 12), whereas they are nwch more 
Irregular below 10 keV (Figure 15). For V = 350 km-s“^, Is 1150 km-s'^, 
and (8.2) Indicates that the turbulent magnetic field spectrum should peak 
near f s ^ 0*04 Hz, or correspondingly down to a wave period of 25 

seconds. Figures 14 and 16 reveal cons1(terab1e p(Mer at about 25-second 
periods when the closed magnetic field bubble sweeps over ISEE-1, -2, and -3, 
and this general behavior persists until the shock. Figure 10 indicated 
that the spectrum peaked between 0.02 and 0.08 Hz, near the shock. We note 
that the 0.01-3 Hz magnetic field variances In Figures 9, 10, 12, and 13 
captured the peak of the wave- spectrum, though they do not contain the wave 
power at periods exceeding 100 seconds produced by protons above 160 keV energy 
Lee (1983) has estimated the rms amplitude of Ion-excited waves at 
the shock front, normalized to that of the upstream Interplanetary field, 
to be about 0.18. Referring to the ISEE-3 measur«nents (which are less 
likely to be contaminated by bow-shock effects) In Figure 11, we find that 
|£B|rms/B was about 0.14-0.16, In excellent agreement with theory. We 


note that the NHD wave amplitudes do not Increase monotonfcally approaching 
the shock, as theory suggests, but appear to fill the entire closed region 
more or less uniformly. Further analysis of the data and extension of the 
magnetic field rms amplitude data to lower frequencies will be needed to 
ascertain vfhether this observation really conflicts with the theory, which 
did not visualize the closed magnetic confinement geometry In which the 
November 11-12, 1978 shock was found. 


8.4 Effects of Energetic Ions Upon Shock Structure 


The November 1978 shock can be used to test ttnK>-f1u1d models 

of shock structure (Drury and Volk, Axford et , 1982) because It 
satisfies the conditions of these theories. We presune that first-order 
Fermi -acceleration theory applies to the strong increase in energetic 
particle fluxes during the last half-hour before shock encounter (see 
Figure 12). This ramp had a scale small compared with that of the magnetic 
bubble, and so we can assimie that the pla»na subshock-energetic particle 
foreshock system was planar. However, significant energetic particle 
fluxes v^re contained by the closed field configuration upstream of the 
extended shock system. In a companion paper (Hada and Kennel, 1983), we 
will extend the above two-fluid shock theories to magnetohydrodynamics and 
compare the predicted energetic particle energy density at the shock with 
observations. We will find that the two-fluid magnetohydrodynamic theory 
overestimates the energetic particle pressure at the shock by a factor 
2-3. Moreover, it Is curious that the slight acceleration of the solar 
wind expected in such theories was not detected In the last half-hour 
before shock encounter. 



8.5 Downstream Properties 


The energetic particle distribution, according to theory, should not 
change at the shock, remain constant In space downstream, maintain the spectrum 
achieved at the shock, and be Isotropic In the solar wind frame. Reinhard ^ 
(1981) and Scholar (1983) have shown that the protons with 10*s of keV 

energy were. In fact, Isotropic In the solar wind frame downstream. Figures 1 
and 12 Indicate that the first three of these requirements were satisfied 
qualitatively by the November 11-12, 1978 shccK (see, also, Reinhard et al . , 
1981). However, the particle Intensity does respond to the passage of the 
rotational discontinuity seven minute: after the shock and maximizes about 
an hour downstream, and the spectrum softens somewhat well downstream of the 
shock. 

Theories of Qy shock structure (Moiseev and Sagdeev, 1963; Kennel and 
Sagdeev, 1967; Auer and Volk, 1973) suggest that thermal anisotropy induced 
by shock expression should create a large-amplitude spectrum of downstream 
transverse electromagnetic waves by the resonant (or non-reso»;ant) firehose 
Instability. A numerical simulation by Kan and Swift (1983) indicates that 
firehose waves will be generated if the Alfven Mach nwnber of the shock exceeds 
about 2.5, a condition satisfied by the Noveni)er 11-12, 1978 shock. Figure 2 
reveals that large-amplitude waves (6By - lOy) existed between the shock and 
the field rotation seven minutes downstream. These appear to have about 10-sec 
periods and a polarization different from that of the waves upstream. Let us 
now estimate the wavelength of these waves. Let us assume that the observed 
frequency Is given by f » KV2/(2w), where V 2 ■ 650 km/sec Is the solar wind 
speed downstream of the shock. With f a 0.1 Hz, we find K a 10“’ cm"^. For 
the measured proton temperature of 1.0x10*® K and a magnetic field of 16 y» 


th« thermal proton Larmor radius Rj Is about >10^ cm, so that KRj a i. Thus, 
we conclude that the periods of the waves observed downstream are consistent 
with the resonant fireshose Instability, whose growth rate peaks near KRj s 1 
(Kennel and Scarf, 1968). 


8.6 Summary 


The observations of the November 11-12, 1978 interplanetary shock 
discussed here by Lee (1983) and by Scholar (1983) confirm the 

conceptual outlines of current theories of quasi-parallel shock structure 
and energetic particle structure. The spatial profiles, energy spectra, 
and angular anisotropies of the shock-associated energetic particles are 
In semi -quantitative accord with first-order Fermi -acceleration theory. 

The frequency and amplitudes of the low-frequency electromagnetic waves 
observed upstream of the shock again agree semi -quantitatively with a quasi- 
11 near description of the resonant cyclotron interaction with energetic protons. 
Large-amplitude turbulence, whose wavelength Is consistent with firehose 
Instability, extended about 50 downstream of the shock. 

It is not clear to what extent the observations of the November 11-12, 

1978 shock are typical. First of all, even moderately strong quasi-parallel 
interplanetary shocks are relatively rare. Moreover, the November 11-12 shock 
appears to have been propagating in a closed magnetic field region when it 
encountered ISEE-1, -2, and -3. This certainly affected escape conditions 
for energetic particles. It does appear that at least the strongest part of 
the energetic proton shock was smaller than the magnetic confinement region, 
and so we used the data from the last half-hour before shock encounter to 
test theory. On the other hand, halo electrons probably had a longer scattering 
mean free path than energetic protons, and their distribution was quite different 
on open and closed field lines. Finally, the electron heat flux upstream of 
Interplanetary shocks may be sensitive to shock parameters as well as geometry. 

Certain aspects of the present study are incomplete. Not much is known 
about few kev "seed" particles upstream of Interplanetary shocks, and although 


we succeeded In detecting them. It requires very high time-resolution 
measurements to document how they emerge from the thermal distribution 
in the shock front. We suggested that these particles might have Interacted 
with the electron heat flux via the Ion cyclotron Instability upstream of 
the November 11-12, 1978 shock, but we have no Indication how common this 
effect might be. Our diagnosis of the turbulent state downstream of the 
shock lacks essential details, such as careful studies of the wave polarizations 
and thermal anisotropies. Such studies appear essential for quasi-parallel 
shocks, which, if they have extended regions of disturbance upstream, should 
also have them downstream. 

Two aspects of the present study are puzzling. The MHD wave amplitudes 
did not build up In the last half-hour before the shock, when both the 
energetic particle fluxes and Ion acoustic wave amplitudes Increased 
significantly. The acceleration of the solar wind expected in the foreshock 
was not detected. 

One question raised by the present study is unsettled; The role of 
the electrostatic and/or electromagnetic electron heat flux instabilities 
in the upstream region. Moreover, we note that our documentation of the 
wave modes present is incomplete. For example, we have not searched for 
the - 1 Hz waves that Sentman et ai. (1983) find are due to energetic electrons 
escaping from the bow shock. 

Clearly, the results of any case study have limited generalities. 

Given that shock structure depends sensitively on the shock parameters, 
and that quasi-parallel shocks are so thick that their properties can depend 
on their global configuration, systematic studies of many shocks are called 
for. We hope the results of this paper can usefully guide such studies. 


9. CONCLUDING REMARKS 


Although theoretical irodels of quas1<>para11e1 shocks are over 20 years 
old, and although the suggestion that energetic particles are significant 
to shock structure Is equally venerable, experimentalists could do little 
with these Ideas until the past five years. The earth's foreshock had a 
complex phenomenology whose disorder had to be reduced before It could be 
fitted Into a theoretical fraiwwork that had seemed peculiarly 111 -adapted 
to bow-shock observations. Now, It Is clear that quasi -parallel shocks 
have such an enormous spatial scale that interplanetary shocks are a better 
experimental arena to test theories of their structure. A coherent viewpoint 
Is now emerging from the experimental and theoretical research of the past 
five years — a viewpoint which this paper is based upon and largely substantiates. 
Only quasi-parallel shocks have the large regions of magnetohydrodynamic 
turbulence upstream and downstream that is an essential Ingredient for first- 
order Fermi-acceleration of energetic particles. This fact is linked to the 
relative ease with which superthermal and energetic particles can stream 
through quasi-parallel shocks. It is clear that such particles generate the 
wave-fields that scatter them, and the outlines of a theory that will eventually 
predict the intensity and spectrum of shock-accelerated particles as a function 
of shock parameters are in view. That said. It Is prudent to add two cautionary 
warnings. First, not much is known about the microstructure of quasi-parallel 
shocks, or even whether they have a microstructure. The current theoretical 
models are based on the Interactions between particles and electromagnetic 
waves with wavelengths equal to or longer than a thermal Ion Larmor radius. 

While It Is conceivable that these physics could account both for the shock 
dissipation and energetic particle dissipation. It Is not proven that It can 


do so. Second, our experience with the earth's bow shock Indicates that the 
shock structure Is strongly parameter-dependent, and the picture of the 
large-scale quasi -parallel shock structure that has energed fr^i the two 
Interplanetary shock studies completed to date might be misleading. It remains 
for future research to confirm, or to temper, our present enthusiasm. 


Figure 1: Global View of the November 11-12, 1978 Shock Event 

The eight panels In this figure show selected ISEE-3 data 
for the period 1800 UT on November 11, 1978, to 2400 UT on 
November 12, 1978. The top two panels present energetic Ion 
measurements; the middle three, selected plasma wave observations; 
and the bottom three, the density, bulk velocity, and magnetic * 
field In the solar wind. Several events are apparent In more 
than one diagnostic. For example, at 2230 UT, the solar wind 
plasma density drops, the 17.8 Hz whistler magnetic amplitudes 
and the 50 kHz thermal continuum amplitudes diminish; the 3.16 
kHz Ion acoustic electric field amplitudes Increase, and the 
35-56 keV proton Intensity reaches a plateau. All these diagnostics 
Indicate that the solar wind entered a new state after 2230 UT. 

The Interplanetary shock, vdilch encountered ISEE-3 at 0028 UT 
on November 12, Is apparent In all the diagnostics. Note that 
the 35-56 keV proton intensity increased by more than an order 
of magnitude in the last half-hour before the shock encounter. 

The 35-56 keV proton Intensity maximized approximately an hour 
downstream of the shock. In rough association with sporadic 
enhancements of the 50 kHz electron plasma wave continuum 
and a burst of ion acoustic waves. Note the sharp drop in 
35-56 keV protor intensity at 1715 UT, which coincided with 
the arrival of the shock driver gas at ISEE-3 (Bame et al . , 

1981). 
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Figure 2 : ISEE-1 and -3 Magnetic Field Measurements at Shock Encounter. 

Shown here are the Y-component (GSM) and the magnitude of the 
magnetic field measured during and after the ISEE-1 and -3 shock 
encounters at 0058:30 UT and 0028:20 UT, respectively, on November 12, 
1978. Note that the time scales are identical, but that the ISEE-1 
and -3 vertical scales differ. The magnetic field magnitude profiles 
are essentially identical at each spacecraft. The shock compresses 
the field from 7.5 to 16. 5y: The shock is followed by a small overshoot 
in the magnetic field; after the overshoot, the field slowly increases 
until the strong rotation, which is apparent in the Y-component, is 
encountered at 0105:20 (ISEE-1) and 0034:10 (ISEE-3). The region 
between the shock and the rotation contains extremely large transverse 
magnetic field fluctuations which hardly affect the field magnitude. 

We will argue that these have wavelengths comparable with the thermal 
ion Larmor radius in the shock-heated plasma and are, therefore, 
consistent with the resonant firehose instability (Moiseev and Sagdeev, 
1963; Kennel and Sagdeev, 1967; Kennel and Scarf, 1968). 
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Figure 3a: Global Geometry of the Interplanetary Shock at the Time of 

the ISEE-3 Encounter: X-Z Plane 

This figure sketches, approximately to scale, the relative positions 
of the terrestrial magnetosphere, ISEE-1, ISEE-2, ISEE-3, and the interplanetary 
shock at 0028 UT on November 12, 1978. Also shown are the magnetic field 
directions measured at ISEE-1 and -3 at that time, together^with the post- 
shock magnetic field measured at ISEE-3. The shock normal n and the 
corresponding shock plane were taken from the Method 3 shock normal solution 
discussed in Section 3.2. The interplanetary field was essentially in 
the X-Z plane, and the shock normal angle 6 bn ^as roughly 40**. The fact 
that the magnetic fields measured at ISEE-1 and -3 were essentially parallel 
indicates that the IMF upstream of the shock was essentially uniform over 
the roughly 200 R£ separating ISEE-1 and -3. A straight-line projection, 
as well as a more sophisticated time-projection (Greenstadt et ^. , 1982), 
of the ISEE-3 field line indicates that it was not connectecTTo the bow shock. 
ISEE-1 and -2, on the other hand, were probably magnetically connected to 
both shocks. 


Figure 3b. Global Geometry of the Interplanetary Shock at the Time of 
the ISEE-3 Encounter: X-Y Plane 

This figure has the same format as Figure 3a. The magnetic field 
rotation apparent in the By-components of Figure 2, followed the shock at 
ISEE-3. If this rotational discontinuity co-moved with the post-shock 
flow (670 km-s“^), it would have been approximately 40 Re downstream of 
the shock and is so sketched. We assumed that both the shock and the 
rotational discontinuity are planar in the absence of information to the 
contrary. Note that the rotation carries the IMF into an anti-gardenhose 
configuration. 
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Figure 4: Electric and Magnetic Fields Before and After ISEE-1 Shock 

Encounter 

The top two panels show the X and Y- components of the interplanetary 
electric field, and the bottom three, the X, Y, and Z-components of 
the interplanetary magnetic field measured at ISEE-1 between 0050 and 
0115 UT on November 12, 1978. J]je third component of the electric field, 

Ez, can be estimated assuming E*B = 0. The measured ExB drift upstream 
of the shock may be shown to be consistent with the measured solar wind 
velocity. The shock at 0058:30 is most apparent in Ey and Bz. Note the 
gentle ramp in E^ that leads the shock by about one minute. The rotational 
discontinuity at 0105:30 is manifested in Ex, Ey. By, and By. Since By 
changes only slightly near 0105:30, we infer that the rotational discontinui 
was propagating in either the + or - X-direction. The large-amplitude 
magnetic turbulence in Figure 2 was restricted largely to the By-component, 
suggesting that it was linearly polarized. This turbulence was confined 
to the region between the shock and the rotational discontinuity. 
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Figure 5: First Magnetic Neutral Sheet Encounter: ISEE-1 and -3 

Both ISEE-1 and -3 encountered a magnetic neutral sheet approximately 
26 minutes after their shock encounters. The top four panels show the 
ISEE-3 magnetic field magnitude and components, and the bottom panel shows 
the ISEE-1 magnitude, for the twelve minutes surrounding the neutral sheet 
encounters. Note that 3x and By change sign at ISEE-3, while Bz remains 
small . 


Figure 6: Next Three Magnetic Neutral Sheet Encounters: ISEE-3 

These three panels show the magnitude of the IMF for twelve minutes 
surrounding the ISEE-3 magnetic neutral sheet encounters at 0114 UT, 0207 UT, 
and 0229 UT on November 12, 1978. The facts that the magnetic field was 
in the X-Z plane (Figures 3a, b) and multiple neutral sheets were encountered 
suggest that ISEE-1, -2, and -3 and the earth were in the general field 
reversal region of the solar wind during the shock event of November 11-12, 
1978. Furthermore, we expect the IMF field topology was complex, with 
alternating regions of closed magnetic loops and "open" field lines. We 
will argue that the shock was in one such loop, which ISEE-3 entered at 
2230 UT. 
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Figure 7: ISEE-3 Electron Model Parameters Upstream of the Shock 

The ISEE-3 solar wind electron data are normally reduced, in survey 
mode, in terms of a model that fits the distribution by a sum of bi-Maxwel 1 ian, 
drifting core and halo distributions. Various moments of these models are 
then displayed as a function of time. The top panel in this figure shows 
the heat flux parallel to the magnetic field, the middle panel shows the 
total and halo component densities, and the bottom panel shows the halo (top 
pair of curves), the core (bottom pair), and the average (middle pair) electron 
temperatures. The thermal anisotropy may be estimated from the width of the 
shaded regions; the top of each region is the parallel temperature, and the 
bottom is the perpendicular. 

Note that the heat flux, which was always away from the sun and the 
oncoming interplanetary shock, fluctuated violently before 1840 UT. The 
heat flux and the densities diminished, and the temperatures increased, 
at 2230 UT and then held constant until shock encounter. The solar wind 
electrons were clearly in a different state after 2230 UT. 
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Figure 8: ISEE-3 Halo Electrons Upstream of the Shock 

The top panel presents the parallel drift speed, relative to the 
solar wind ions, of the ISEE-3 model halo electron distribution. The 
middle and bottom panels present the partial densities and first moment 
drift velocities of the away (V|| > 0) and toward (Vn < 0) semi-distributions 
of halo electrons. The net halo drift (top panel) diminished near 2230 UT, 
not because the number of electrons streaming away from the shock declined, 
but because they were almost compensated by a new component of electrons 
streaming toward the shock. Since the ISEE-3 field line was not connected 
to the bow shock, the new "toward" halo electrons cannot have come from the 
bow shock. We believe that they, too, were generated at the interplanetary 
shock and streamed the long way around a closed magnetic loop to ISEE-3. 

This loop had a spatial scale of 800 R£ or more. 





Kigjre 9: ISEE-3 Entry into the Closed Magnetic Loop 

The increase in 3.16 kHz ion acoustic electric field amplitude at 
2230 UT in Figure 1 evidently was associated with ISEc-3's passage ■‘nto 
the closed magnetic loop. The bottom panel in this figure shows 1-min 
average 3.16 icHz electric field amplitudes on an expanded time scale. The 
averages clearly increase around 2230 UT. The top panel indicates that 
the 35-56 keV proton flux increased about a factor 3 upon entering the 
loop. The middle panel shows the first parallel anisotropy moment of the 
energetic protons. The protons become essentially isotropic upon entering 
the loop. Figure 1 indicates that there was sporadic activity in the ion 
acoustic amplitudes prior to 2230 UT; however, the intensities remained 
high after 2230 UT and increased toward shock encounter. 
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Figure 10: Upstream Ion Acoustic and MHO Waves: ISEE-3 

The top panel shows the 3.16 kHz ion acoustic wave amplitudes, and 
the bottom two panels show the variances in the .01-3 Hz interplanetary 
By and Bz field components for the period 2100 UT on November 11 to 0100 UT 
on November 12, 1978. The times of mognetic bubble entry, and the encounters 
with the shock, field rotation, and the first neutral sheet are Indicated 
by vertical lines. Bursts of ion acoustic and MHO wave activity accompany 
one another before bubble entry, and both ion and MHO waves are continuously 
enhanced within the bubble. Note that the 1-min average ion wave amplitudes 
increase significantly in the last half-hour before the shock, during which 
time the 35-56 keV proton intensity also increased, whereas the .01-3 Hz 
magnetic field amplitudes remained more or less constant. 
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Figure 11: Frequency, Polarization, and Amplitude Characteristics of MHD 

Waves Measured Upstream of ISEE-3 Shock Encounter 

Shown here are the wave normal angle SkBi the ellipticity, the percent 
polarization, and the total rms and compressional amplitudes (normalized 
to the magnitude of the IMF) plotted against the frequency measured in the 
spacecraft frame. The spectral decomposition between .005 and 0'.32 Hz for 
the period 0012-0028 UT on November 12, 1978 is plotted as a solid line. 

An indication of the spectral properties above 0.32 Hz is given by the 
dotted curves, which show data for the period 0025-0027:30 UT which ended 
approximately 40 seconds before shock encounter. 

The total rms power was enhanced between 0.01 and 0.16 Hz; the proton 
cyclotron frequency was roughly 0.12 Hz. The absence of significant power 
near 1 Hz indicates that the wavelengths were longer than the thermal ion 
Larmor radius. Note that the waves were relatively weakly, elliptically 
polarized and propagated more or less parallel to the IMF (6 kb • 20°). The 
peak-normalized rms amplitude (0.16) will be compared with theoretical expectati 
in Chapter 8. 



Figure 12: Energetic Protons and MHD Turbulence Upstream of the ISEE-3 

Shock Encounter 

The top panel shows the energetic proton fluxes in eight different 

energy channels; the middle panel, rough spectral indices formed from the 

ratios of the fluxes in adjacent energy channels; and the bottom panel, the 
variance in the Y-component of the IMF for the period 2100 UT to 0100 UT 
on November 11-12, 1978. The channel energies are as follows: (1) 35-56 keV; 

(2) 50-91 keV; (3) 91-147 keV; (4) 147-237 keV; (5) 237-383 keV; (6) 383-619 keV; 
(7) 619-1000 keV; and (8) 1000-1600 keV. The scales appropriate to the energy 
channels alternate: Channel 1 to the topmost left, channel 2 to the topmost 

right, and so on. The labels 1/2, 3/4, etc. in the middle panel refer to the 
energy spectral index inferred from the ratios of the fluxes measured in the 
first and second channels, third and fourth, and so on. 

Note that the fluxes in channel 1 increased before bubble entry, 

plateaued, and then diminished near 2330 UT. The protons between 35 keV 

and 237 keV increased in intensity after 0000 UT and then held roughly constant 
after the shock. Protons at all measured energies responded to the rotational 
discontinuity. The fluxes of protons above 619 keV (channels 7 and 8) rose 
perceptibly as the shock approached, but did not change as the magnetic bubble 
passed over ISEE-3. The low-energy spectral indices diminished as the shock 
approached, and the 1/2 index even became negative, indicating the presence 
of high-energy particles diffusing to ISEE-3 from the oncoming shock. 











Figure 13: Plasma Wave Electric Fields Observed Upstream of the Interplanetary 

Shock and Bow Shock: ISEE-1 

Presented here are swept frequency receiver (SFR) plasma wave electric 
field measurements made between 1800 UT and 0200 UT on November 11-12, 1978. 
ISEE-3 exited the bow shock at 1840 UT and immediately began to detect 
intense 40-50 kHz electron plasma waves which persisted, with a 'decrease 
in frequency consistent with the general decrease in density measured at 
ISEE-3 (Figure 1), until the interplanetary shock at 0058 UT. Such plasma 
waves were not detected at ISEE-3, which was not magnetically connected 
to the bow shock. The fact that the plasma waves at ISEE-1 were detected 
from bow-shock exit until interplanetary shock encounter suggests that 
ISEE-1 was connected to the bow shock, and probably also to the interplanetary 
shock, as Figures 3a, b suggest. 





Figure 14: Comparison of Upstream MHD Waves Measured at ISEE-1 and -3 

The top two panels show the RMS variances in the magnitude, Bj, and 
the Y-component, By. of the interplanetary field at ISEE-3 for a 3i-hour 
period surrounding shock encounter. The bottom two panels show the same 
quantities measured at ISEE-1. We estimate that ISEE-3 entered a closed 
magnetic loop at 2228 UT on November 11; if the loop convected with the 
solar wind speed between ISEE-3 and -1, ISEE-1 would have entered it at 
about 2330 UT, as indicated. The variances at ISEE-1 do increase at that 
time. Note that the ISEE-1 magnitude variances are considerably larger 
than those at ISEE-3. This difference may be due in part to the fact that 
ISEE-1 may have detected upstream turbulence from both the interplanetary 
and bow shock. 





Figure 15: Bubble Entry and Shock Encounter: ISEE-1 

The top two panels show 1 kHz and 1.78 kHz ion acoustic electric field 
amplitudes for the period 2200 UT - 0200 UT on November 11-12, 1978. The 
1-min average amplitudes are indicated by solid shading, and the peak amplitudes 
detected during the averaging interval are connected by a thin line. Large 
peak-to-average amplitude ratios characterize solar wind ion acoustic 
turbulence. The middle panel shows the variances in the .01-3 Hz Z-component 
of the IMF, and the bottom two panels show the differential fluxes of 2 and 
6 keV protons. Upstream of the bubble entry near 2330 UT, the activity in 
ion acoustic waves, MHO waves, and few keV protons tends to be correlated 
in time. Bubble entry is marked by a striking increase in the 6 keV proton 
differential flux. 
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Figure 16: Character of Upstream MHD Waves at ISEE-3 and ISEE-1 Magnetic 

Bubble Entry 

Although MHD wave activity occurred at lSEE-1 and -3 prior to bubble 
entry, it was continuously present afterward. This figure compares the 
MHD wave activity in all three vector components and in the magnitude of 
the magnetic field for a 12-minute period surrounding ISEE-3 bubble entry 
at 2228 UT, and the ISEE-1 entry at 2328 UT. ISEE-3 data is at the top 
of each panel, and ISEE-1 data is at the bottom. The ISEE-3 magnitude scale 
is at the left of each panel, and the IStE-l scale is at the right. The 
ISEE-3 time scale is at the top of the figure, and the ISEE-1 scale is at 
the bottom. It is clear that a significant enhancement of wave activity 
occurred at bubble entry and that the waves had the same general periodicity 
and character at both spacecraft. In particular, the magnitude and Y-component 
was noticeably less variable than were the Y and 7.-components of the 
interplanetary field. The fact that large-amplitude oscillations commence 
within a few wave periods suggests that the waves were already present in 
the bubble before it swept over the spacecraft. 
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figure 17: Modulated Few keV Protons at Bubble Entry: ISEE-1 and -2 

This figure compares the 6 keV proton fluxes and the Z-component 
of the IMF measured at lSEE-1 and -2 for a 12-minute period surrounding 
bubble entry. Extremrly large modulations of the proton fluxes accompany 
the MHD wave activity; the peak fluxes occur when B? Is closest to zero, as 
is indicated by the vertical lines. Since the 6 keV proton detector looks 
essentially along the Z-direction, the peaks correspond to 90* pitch angle 
protons. The 6 keV protons appear to be extremely anisotropic. 

ISEE-2 was about 4000 km upstream of lSEE-1, corresponding to a 12-second 
delay at the solar wind speed. The ISEE-2 »'ata has been offset relative 
to the ISEE-1 data by 12 seconds. The rough correspondence of the proton 
intensity peaks and magnetic field oscillations suggest that protons and 
wave' were convected from ISEE-2 to ISEE-1. 
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Figure 18: Electron Drifts and Ion Acoustic Waves: ISEE-3 

The large electron heat fluxes in Figures 6 and 7 would lead to a strong 
parallel current unless the core electrons drifted relative to the ions in a 
direction opposite to the halo electron drift. The top panels compare the 
electron current carried by the halo electron model distribution with the 
current carried by the core electrons. The extent to which the -two profiles 
are mirror symmetric indicates the degree of current neutrality. Note that 
the core electrons drifted with a oarallel speed of about 200 km-s"^ after 
bubble entry at 2228 UT. The bottom panel shows 5-min average (bar graph) 
and peak (dots) 1.78 kHz ion accoustic electric field amplitudes for the 
same time period. The origin of solar wind ion acoustic turbulence has been 
puzzling, because it is generally believed that the electron-to-ion temperature 
is too low for instability during such events. Indeed, the 200 km-s"^ measured 
drift was below ion acoustic threhold for this event; it is possible that the 
ion wave activity was thermal fluctuations enhanced by the large electron 
heat flux. On the other hand, the measured core dr’Ft was near or at threshold 
for the ion cyclotron heat flux instcjility. 



